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Abstract: Dressed Gluon Exponentiation (DGE) is used to calculate the thrust and the 
heavy-jet mass distributions in e~^e~ annihilation in the two-jet region. We perform a 
detailed analysis of power corrections, taking care of the effect of hadron masses on the 
measured observables. In DGE the Sudakov exponent is calculated in a renormalization- 
scale invariant way using renormalon resummation. Neglecting the correlation between 
the hemispheres in the two-jet region, we express the thrust and the heavy-jet mass 
distributions in terms of the single-jet mass distribution. This leads to a simple descrip- 
tion of the hadronization corrections to both distributions in terms of a single shape 
function, whose general properties are deduced from renormalon ambiguities. Matching 
the resummed result with the available next-to-leading order calculation, we get a good 
description of the thrust distribution in a wide range, whereas the description of the 
heavy-jet mass distribution, which is more sensitive to the approximation of the phase 
space, is restricted to the range pn < 1/6. This significantly limits the possibility to 
determine as from this observable. However, fixing ag by the thrust analysis, we show 
that the power corrections for the two observables are in good agreement. 
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1 Introduction 



Hadronic jets from the fragmentation of high-energy partons is one of the most spectac- 
ular phenomena in particle physics. It also provides evidence for the partonic nature of 
hadrons: observables such as event-shape distributions can be calculated perturbatively 
in terms of quarks and gluons and compared with hadronic final state measurements 0]. 
At the same time, these observables are directly sensitive to the hadronization process, 
providing an opportunity to learn more about the confinement mechanism from experi- 
ment. 

The LEP experiments and their high quality data have prompted much theoretical 
progress in the understanding of QCD radiation and in the development of appropriate 
tools. This includes Sudakov resummation [0-0, renormalon resummation and 
parametrization of power corrections |7[]- [p^ . 



Our goal is to extract from QCD quantitative predictions for event-shape distributions 
in the two-jet region. The precise LEP data provide a strong motivation for such an effort: 
it can be used for precision measurements of a^, but most importantly to test and improve 
our understanding of the QCD dynamics. An essential requirement is to understand how 
perturbative QCD can provide a reliable starting point for the calculation of event-shape 
distributions. Only then can one address the question of non-perturbative corrections. 

The region where perturbation theory applies is limited from both ends. The multijet 
region is in principle accessible by going to ever higher orders. In practice, however, these 
calculations soon become difficult. On the other hand, the cross section is large only in 
the two-jet region, where it is determined by soft and collinear radiation. Close enough 
to the two-jet limit, the relevant scale for gluon emission becomes of order of the QCD 
scale A, and the problem becomes essentially non-perturbative. 

In the two-jet region one must take into account the two immediate consequences of 
the restricted phase space, namely multiple emission and low scales. Even a qualitative 
description of the distribution da/dy {Q, y), where Q is the centre-of-mass energy and y is 
an event-shape variable that vanishes in the two-jet limit, requires a//-order resummation 
of Sudakov logs L = Inl/y. Such resummation guarantees the vanishing of the cross 
section at y = 0, whereas any fixed-order calculation diverges in this limit. 

Sudakov logs can be resummed thanks to the factorization property of QCD matrix 
elements. The minimal resummation required for quantitative predictions is one where 
the Sudakov exponent is calculated to next-to-leading logarithmic (NLL) accuracy 
In NLL resummation the exponent includes terms of the form a'^L^'^^ and a^L" for any n, 
while terms of the form a^U^ with m < n are neglected. The NLL exponentiation 
formula has become the standard tool for data analysis in recent years. Superficially 
this resummation applies for arbitrarily small y. However, due to the running of the 
coupling it actually applies only for yQ 3> A. In particular, when yQ ~ A not only 
does the perturbative approximation break down completely, but also power corrections 
of arbitrary power ~ (A/(yQ))" become important. 

Since the distribution peak is located at yQ > A, Sudakov resummation with a fixed 
logarithmic accuracy, such as NLL, is insufficient there. The low scales involved manifest 
themselves in renormalons: large perturbative coefficients that increase factorially at large 
orders. These large coefficients enhance the subleading logs a"L'" with m < n, disrupting 
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the dominance of the leading logs. Estimates |^ of the contribution of subleading logs, 
which are neglected in the standard NLL resummation, to the thrust (T) distribution at 
(1 - T)Q > A is of the order of 20% at Q = Mz. 

In Dressed Gluon Exponentiation (DGE) ^ the Sudakov exponent is calculated, 
starting with an off-shell gluon emission. Following the spirit of BLM @-[^] and [p^ , |^, 
the actual gluon virtuality is used as the scale of the running coupling. Through the 
integral over the coupling DGE resums the renormalon contribution (in the large /?o 
limit) to the Sudakov exponent, and thus it takes into account the factorially enhanced 
subleading logs to any logarithmic accuracy. As a consequence, the calculated exponent 
is renormalization-scale invariant. 

The calculation of the exponent by DGE makes the simplifying assumption that glu- 
ons are emitted independently. While reproducing the exact result to NLL accuracy [0, 
this simplification amounts to approximations in both the phase space and the matrix 
elements beyond this accuracy. Nevertheless, since the two sources of large coefficients, 
Sudakov logs and renormalons, are accounted for, we consider DGE a reliable start- 
ing point for the analysis of the distribution in the range yQ > A. Eventually, the 
dominance of the corrections we resum should be tested when the complete next-to-next- 
to-leading logarithmic (NNLL) and next-to-next-to-leading order (NNLO) corrections 
becomes available. 

Physically it is clear that non-perturbative corrections due to soft gluon emission 
are important in the two-jet region. Technically, the DGE calculation gives a clear 
indication for the presence of power-suppressed corrections: the resummation of infrared 
renormalons is ambiguous at power accuracy. Consistency of the full cross section implies 



the existence of genuine power corrections of the same form. As in ||TT|, IT2[ and |]T9| 



we will assume that these power corrections dominate over power corrections that have no 
trace in the perturbative calculation. This has far-reaching implications concerning the 
functional form of the non-perturbative correction. The leading corrections can be written 
as a convolution of the perturbative result with a shape function of a single variable yQ, 
as suggested by Korchemsky and Sterman |jl3|-||T^, 0. The simplest approximation to 
the shape function amounts to a shift of the perturbative distribution, in accordance 
with |T^ , [23(1 . Based on the renormalon ambiguity (in the large /5o limit) only power 
corrections corresponding to odd central moments of the shape function are present. 
This gives, on the one hand, an explanation for the success of shift-based fits and, on the 
other hand, an opportunity to test the renormalon dominance assumption by fitting the 
data and extracting higher moments. 

The DGE method has been the basis of our analysis of the thrust distribution in §]. 
The comparison with experimental data was encouraging in several respects: first, good 
fits were obtained in a wide range of energies and thrust values. These fits were quite 
stable under a variation of different parameters, such as the cuts on the fitting range 
in thrust. In addition, the value of and the non-perturbative parameter were consis- 
tent between the shape-function- and the shift-based fits. Moreover, these values were 
also consistent with those obtained by fitting the average thrust using renormalon re- 
summation [0]. Such consistency, which is a crucial test of whether the approximation 
used is satisfactory, is not achieved ^1|] in standard NLL fits where renormalons are not 
resummed. 
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One aspect that certainly deserves attention is the fact that the best fit value of a 
obtained in 0, ^ (Mz) = 0.110, is significantly lower than the world average p2 
value af^ (Mz) = 0.1172 ± 0.0020. The experimental uncertainty is very small and the 
sources of theoretical uncertainty investigated in hardly sum up to 5%. The latter 
uncertainty is dominated by non-logarithmic corrections at NNLO, which should become 
available in the foreseeable future. 

It is of great interest to extend the application of DGE to more event-shape distri- 
butions. The main motivation is to test and improve our understanding of various per- 
turbative and non-perturbative questions. In particular, it is interesting to investigate 
the relations between non-perturbative corrections to different observables. The idea of 
a "universal" (observable-independent) origin of the leading power corrections to various 
event-shape observables has been the subject of much theoretical and phenomenological 
activity in the past few years. Universal power corrections were suggested both within 
the infrared finite coupling approach |]l9l-^T|, and in the shape-function approach 



In either case, the observable independence of the power corrections involves additional 
assumptions. 

Probably the most relevant pair of variables to be compared in this respect are the 
thrust and the heavy-jet mass. For completeness we recall their definitions here, together 
with some other, closely related, jet mass variables: 



Thrust T 



Heavy-jet mass 
Light-jet mass 



Hi \Vi ■ 
max — - — 

Ei \Vi\ 



Single-jet mass pm = 



pH = max(pi,p2) 
PL = min(pi,p2)- 



(1) 

(2) 

(3) 
(4) 



Here the index i runs over all the particles and m = 1 , 2 represents the two hemispheres 
B.jn defined by the thrust axis ut. 

There are several useful relations between these observables, which follows from their 
definition: the single-jet mass is the average of the heavy and the light. 



1 da 
a dp 



1 da 1 da 

+ 



a dpH (y dpL 



(5) 



In the two-jet region t = 1 — T equals, up to corrections suppressed by t, the sum of the 
two jet masses 

pH + pL ■ (6) 

To leading order in as this relation is exact and, since the light-jet mass vanishes at this 
order, it follows that 

1 da 1 da Ida /„s 
— r = -^ = 2-^ {0{as)). (7) 
a dt a dpH <J dp 

In the approximation of independent emission used here, both the thrust and heavy-jet 
mass distributions can be expressed in terms of the single-jet mass distribution . This 
imphes that a single shape function can be used to parametrize the non-perturbative 
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corrections in both cases. On the other hand, fits based on the NLL resummation formula 
and a shift for the two distributions yield different values of a<j, along with a different 
magnitude for the power correction (the shift) |31[]. While a fair comparison of the 
power corrections must assume a common value of as-, such a fit does not improve the 
agreement between the power correction parameters. We will see below (Fig. ^ that 



fitting the heavy-jet mass distribution by shifting the NLL distribution, as in |31], but 



fixing the value of based on the thrust fit, the power correction of the heavy-jet mass 
is much smallei^ (~ 50%) than expected based on the thrust fit. 

Since the relation between the power corrections to the heavy-jet mass and those 
to the thrust followsf] from the definition of these variables, the discrepancy between 
the power corrections must be interpreted as a sign that the approximations made in 
the calculation are too crude. Based on |^, 0|, it is quite clear that one reason for this 
discrepancy is the fact that renormalon contributions were neglected. We will verify this 
assertion in what follows. 

In , the relation between the non-perturbative corrections to the two distributions 
was explained in terms of a strong correlation between the two hemispheres due to non- 
perturbative soft gluon emission at large angles. It was shown that introducing a positive 
correlation through a two-hemisphere shape function, the heavy-jet mass and the thrust 
(as well as the C-parameter) can be fitted. However, such a correlation, if important, 
should appear first at the perturbative level. Since the resummation formulae implicitly 
assume independent emission, and thus no correlation, we would rather not introduce 
such a correlation on the non-perturbative level. We shall assume here that in the two- 
jet region the jet masses are uncorrelated also on the non-perturbative level, and confront 
this assumption with the data. 

Another important issue in the comparison of power corrections to the thrust and 
the heavy-jet mass is the dependence of the distribution on the hadron masses |2^. The 
perturbative calculation is based on massless partons, whereas the measured observable 
is based on massive hadrons. It is clear that the spectrum of hadrons has some effect 
on the distribution. In general at high enough Q and y, the details of the spectrum do 
not matter. Instead, some generic characteristics of hadronization dominate, and a small 
set of non-perturbative parameters controlling power corrections is sufficient. However, 
the comparison of the thrust and the heavy-jet mass is problematic since the thrust is 
(conventionally) defined in terms of the three-momenta whereas the heavy-jet mass in 
terms of the four-momenta. This implies that the two observables will exhibit different 
sensitivity to the masses of the measured particles. Ways to overcome this difficulty by 
appropriate modifications of the definitions (which coincide with the standard ones in the 
massless limit) were suggested in [^. We will see that choosing a common modification 
of the definition (hadron mass scheme) for the two observables is indeed crucial for a 
meaningful comparison of the power corrections. 

The outline of this paper is as follows. We start in section 2 by recalling the pertur- 
bative calculation of the thrust and the heavy jet mass distribution using DGE. We then 
discuss the limitations of the approximation at hand, and, finally, we present our assump- 



tThe NLL fits in yield a lower value of for the heavy-jet mass and a higher value of the power 
correction (shift) compared to the thrust values. 

■I- One does not rely here on the assumption of universality of an infrared finite coupling. 
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tions concerning hadronization corrections. Section 3 summarizes the data analysis. We 
begin by performing separate fits to the thrust and the heavy-jet mass distributions and 
then study the consistency between the two. Section 4 contains our conclusions. 



2 Assumptions and calculation 

In this section we present the calculation of the distributions, emphasizing the main 
assumptions involved. In section 2.1 we calculate the jet mass distribution from a single 
dressed gluon (SDG). In section 2.2, we exponentiate the SDG cross section arriving 
at the DGE formula for the single-jet mass distribution. This distribution is then used 
to express the thrust and the heavy-jet mass distributions. In section 2.3 we discuss 
the limitations of the considered approximation. Finally, in section 2.4 we present our 
approach concerning the parametrization of hadronization corrections and the way we 
deal with hadron mass effects. 



2.1 Single Dressed Gluon cross section 

We begin by calculating the differential cross section for a single gluon emission. The 
calculation is performed with an off-shell gluon fc^ = > in order to provide the 
basis for renormalon resummation ^7\, ^ 0, |^, ^. Contrary to @, ^, which used the 
full matrix element and identified the log-enhanced terms prior to the integration over 



the gluon virtuality, we will follow the approximation suggested in p5|], which allows to 
isolate the relevant terms right at the beginning. 

By performing the calculation with an off-shell gluon we will effectively be integrating 
inclusively over the decay products of the gluon. Based on the definitions of the variables 
(E) ^ (i) such a calculation yields the exact result when the partons that the gluon 
decays into are confined to a single hemisphere, but not so when these partons end up 
in different hemispheres. Therefore, our result should be considered as an approximation 



even in the large (3o (large Nj) limit ||3^. There are strong indications that this inclusive 
approximation is good in the case of the thrust [^, 1^, ||, We will return to discuss 
this point in the section 2.3. 

Let us denote the final momenta of the two quarks by p and p, where = = 0. 
Sudakov logs are associated with the singularity of the quark propagator prior to the 
emission, + + kY- Concentrating in the two-jet region, we will neglect non- 

logarithmic terms in the cross section. This is equivalent to neglecting non-singular 



terms in the matrix element It is convenient to use the light-cone axial gauge p5 
where, in the approximation considered, the gluon effectively decouples from one of the 
quarks. 

Assuming that the gluon is in the hemisphere of the quark p, it is p that sets the thrust 
axis (we will calculate the cross section only in this half of phase space and account for 
the other half by an overall factor of 2). We define the light-cone coordinates such that 
p = (p+,0,0) and p = (0,p_,0), and use the light-cone gauge where = 0. In this 
gauge the gluon coupling to the p quark is suppressed^], so it is enough to consider the 



^In the on-shell case, = 0, this couphng is strictly zero: yl+ J_ — 0. 
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diagram where the gluon is emitted from p. The gluon propagator is given by 



V - -9i,v H -j- • (Sj 

The cross section in the centre-of-mass frame is 

The leptonic tensor L'"^ = |Tr {ii~i''i2l'^} can be replaced by \Q'^g''" ■ To calculate the 
hadronic tensor Qp^RP" = Z^spins'^-^^ W6 write the amplitude 

= (^:^^^'^(^')^''(^ + ^)Yu^'\p) (9) 

where t** is a colour matrix and e^* is the gluon polarization vector. The squared matrix 
element summed over the gluon and quark spins is then given by 

^ -MA^^ = d,,TT {(^ + + m'^r} (10) 

spins \ ^ Jr ) 

where '}2s'^^^\'P)^^^\'P) "was replaced by ]/, and X^a ^y*^^ by d^^. Introducing the Sudakov 

decomposition oik = fcj,) = I3p+ap+kj^, and defining the following dimensionless 
parameters 

j3 = kj^/pj^ = 2kp/2pp 

a = k_/p^ = 2kp/2pp (11) 

7 = kj_/2pp 

A = /2pp = {2k^k^ — k\)/2pp = ajS — 7, 

the propagator is proportional to 1/(q; + A)^. As observed in to calculate the log- 
enhanced cross section, it is sufficient to take into account the singular terms in the 
matrix element, corresponding to the situation where both a and A are small (with no 
specific hierarchy between them). From the condition aj3 = 7 + A it follows that the 
transverse momentum fraction 7 is also small (at least with respect to /3) while j3 can be 
either large (the coUinear limit) or small (the soft limit). In these variables 



spins 

The cross section is therefore]]]. 



■A 1 ;3^ + 2/J + 2 



(a + A)2' a + A (3 



(12) 



CpO-s [ dad(3 dX5{\ — e{\ + 13)) 



TT 



a + A)2^ ' a + X (3 



(13) 



^An overall factor of 2 accounting for the other half of phase space is included. 
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where ctq = iVclZ^g e^) 47rQ;^/(3(5^) and the 6 function guarantees that the integration is 
performed for a fixed e = vn^jQ'^ and . The integrand has the interpretation of a 



generahzed, off-shell gluon splitting function p5 ]. 

We are now ready to calculate the differential cross section for the jet mass distri- 
bution. The hemisphere of the quark p to which the gluon is emitted is the heavy one, 
= {p + kY/Q'^, whereas the other has zero mass, m\ = 0. The relation between 
Ph = m'jj/Q^ and the Sudakov variables (|lT]) is simply 



PH 



(p + kf 2pk + (a + A) 2pp 



a + X 



a + X 



1 + a + p + X l + (3 



(14) 



Note that also t = l-T= {mjj + m\)/Q'^ = pn- 

Changing variables in (|13D we get the following result for the single off-shell gluon 
differential cross section 



1 da 



CpOts 



a dp 



H 



TC 



d(3 



-— + 

Ph 



(15) 



Next, in the integration over /3, it is the limit of vanishing transverse momentum 
= X that generates the log-enhanced terms [^. We therefore integrate over P and 
substitute the limit P = e/{pH — e). The result is 



1 da 



Cpas 



a dp 



H 



TX 



2 PH 3 1 e 1 

- In + — + 77 — 

Ph e 2pH Ph ^Ph. 



(16) 



Using the gluon bremsstrahlung effective charge p3[ — > ds is sufficient to guar- 
antee that the singular terms in the splitting functions are exact. This way the renormal- 
ization scheme is not arbitrary: there is a unique effective charge that generates correctly 
subleading terms. Provided one also uses the two-loop renormalization group equation, 
the calculation of the thrust and the heavy-jet mass distributions will be exact ^ to 
NLL accuracy. 

We now replace by an integral over the time- like discontinuity of the coupling p2, 
27| , |28| . In this way we take into account the relevant physical scale for the gluon emission 
/c^ = eQ^. Defining A{k'^) = /3oQ;s(fc^)/7r, we will use here the scheme invariant Borel 
representation p5|-p7|. 



A{k'^) = j^ rfw exp (^-MlnA;VA^) ^b(m), 



(17) 



where Ab{u) = 1 for the one- loop coupling, provided A is defined appropriately. The 
integrated time-hke discontinuity is [0, |2! 



A^g{Q^) = J du exp \^-u\nQ'^/A 



^'^^^ Arm. 



10 ^ ^ TTU 

After integration by parts the single dressed gluon cross section is 



(18) 



1 da 
a dpH 



{Q\ph) 



Cf fpH de 



SDG 



Ph Ph Ph 



(19) 
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Here, the upper integration limit is determined from the condition that a is positive. 
Along the zero transverse momentum boundary of phase space, the limit is a = i^pu — 
e)(l + /9) = 0. In this situation both k\ and vanish, so the gluon is exactly coUinear. 
The lower integration limit is determined from the requirement that /? > a (otherwise the 
gluon belongs to the p hemisphere). Along the zero transverse momentum boundary a/3 = 
A, this condition translates into /5 > \f\. It follows that the lower integration limit is 
e = p'jj. This limit corresponds to large-angle, soft emission. 

To obtain the single-jet mass distribution we have to average over the heavy and the 
light hemispheres. But since the light-jet mass cross section starts at two gluon emission 
(order a^) this amounts, according to (|^, to dividing the heavy jet mass cross section 
by 2, 



a dp 



P de ^ 



p^ p-^ 



(20) 



Substituting (|T8D into (|20|) and performing the integral over e yields the Borel represen- 
tation of the logarithmically enhanced SDG cross section. 



Ida 

a dp 



{Q\p) 



SDG 



Wo 



r duB,,^{u,p) exp (-ulnQ'/A') ^-^^Ab{u) (21) 

Jo ^ ^ ITU 



with the following Borel function 



Bsog{u,p) 



1 


"2 






— exp 


(2„ lui) - ( 


p 


u 





2 1 

u 1 — u 



+ 



u 



exp 



u In 



(22) 



This result was discussed in detail in [|]. Here we just recall that the sensitivity to soft 
gluon emission does not appear here through Borel singularities (which cancel out by the 
sin(7r'u) factor) but rather as convergence constraints on the Borel integration at m — oo. 



2.2 Exponentiation 

Assuming independent emission and additive contribution by each emission to the jet 
mass: p = Ylk=iPk, the multiple gluon phase space factorizes in Laplace space 0. As 
in 0, the DGE cross section is obtained by exponentiating the SDG cross section, given 
by ( pip and (P^, under the Laplace transform. The Laplace transform J(i/, Q^) of the 
single-jet mass distribution is then given by 

lnJ{u,Q')= f'-^{Q\p) 
Jo a dp 



-vp 



1) dp, 



(23) 



where the Laplace weight e is associated with real emissions and —1 with virtual 
corrections. The corresponding Borel representation ^ is 

2(3oJo 

with 

1 



duBiy^u) exp {^—uhiQ'^/k' 



Ab{u) 

iru 



(24) 



B u,u) 



1 dp 
p 

2u In i/p 



^ 2m In i 

— e P 

U 



-2u) + 



u 

1 

2m 



1 



u 
2 

u 



u 



1 



u 



e p 



1 

2-u 



(e-- - 1) 



(26) 



I In V 



Ti-u) + 



u. 



8 



where non-logarithmic terms have been neglected. Contrary to the SDG result (P^), the 
DGE result has a rich renormalon structure, notably, it has singularities at half-integer 
values of u. These Borel singularities indicate power corrections. Note that the first 
term in p5|), which has a double argument 2u, corresponds to large-angle soft emission, 
whereas the second corresponds to collinear emission. 

The single-jet mass distribution can be obtained by the inverse Laplace transform 

of 



1 da 
a dp 



(P, Q') = ^l ^ exp {up + In J(z/, Q2)| ^ ±R^{Q\ p) ^ j(p, Q 
dp Jc zmiy ^ ^ dp 



d 

dp 



(26) 



where the contour C goes from —too to +ioo to the right of all the singularities of the 
integrand. Using this distribution, the thrust t = pi + p2 distribution is 



1 da 
a dt 



{t,Q^) = I dpidp2j{puQ^)j{p2,Q^)Sipi+ p2-t) 

P20(P2 — Pi) distribution is 



(27) 



d 

dt Jc 



whereas the heavy-jet mass pn = Pi0(pi 
1 da 



P2) 



adpH 



J dpidp2jipi,Q'^)jip2,Q'^) [5(pi -p//)e(pi -P2) +5(p2 -p//)0(p2 -Pi)] 

J dpidp2j{pl, Q^)jip2, Q^) Q{pH - Pl)Q{pH - P2) 



dpH 

d 



dp 



H 



d 

dpH 



(28) 



Equations (pT]) and (^) are our final expressions (prior to matching with the NLO 
result) for the thrust and the heavy-jet mass distributions. Using In J(i/, Q^) in (^) 
with the two-loop gluon bremsstrahlung coupling, as we do below, these results are exact 
to NLL accuracy. However, they differ from ^ by large subleading logs which emerge 
from the integral over the running coupling. In order to understand how significant these 
additional terms are, it is useful to rewrite In J(z/, Q^) as an asymptotic series in the 
coupling. 



In J(z/,Q^ 



2^0,-=! 



(29) 



where ^ = AiQ"^) Inv. Truncation of this series amounts to a fixed logarithmic accuracy. 
The first two functions {k = 1, 2) coincide with the known leading and next-to-leading 
log functions, whereas fk for A; > 2 are given by 

/3(0 = 0.804/(1 -0 -2.32/(1 -20 
/4(0 = 0.779/(1 - 0' - 2.68/(1 - 20' 
/5(0= 2.32/(1 -0^-5.00/(1 -20^ 
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where, for simplicity, we used a one-loop running coupling. At higher orders, these 
functions behave as 



m - [(Y3^-2-5^Y^^J- (30) 

The factorially increasing coefficients of fk{0 enhanced singularity at ^ = 1/2 

imply that a fixed logarithmic accuracy calculation holds only in a very restricted range 
in the A{Q'^) and L = Inz/ ~ Inl/p parameter space. Note that the contributions of 
two consequent terms in ( ]29|) are of the same order when A(Q^) ^ (1 — 2^)/{k — 2). At 
this order the perturbative expansion is exhausted (the minimal term is reached) and 
the correction is effectively power-like ~ A/{Qp). We assume below that the order of 
the minimal term is higher than NLL, so including higher orders in (^) improves the 
evaluation of In J(i/, Q^). Figures 2 and 3 in show that near the distribution peak, and 
anywhere above it, p > p(peak), this is indeed the case. The fixed logarithmic accuracy 
calculation is justified if ( pOf ) converges well. For a fixed C, (any ^ < 1/2) this holds if 
the coupling is very small ^(Q^) ^ 1 — 2^. However, for a given coupling, there is good 
convergence only for ^ -C 1/2, i.e. 



p > exp 



A 



2A{Q^) Q- 



This condition is not realized in the peak region at any relevant centre-of-mass energy. 
Consequently a fixed logarithmic accuracy is insufficient. The analysis of showed, 
that in the case of the thrust, the contributions of NNLL and further subleading logs to 
In J(z/, Q^) at Q = Mz lead to a correction of ~ 20% for p > p(peak). Only when all the 
terms (up to the minimal term) in the series (^Uj) are kept can one safely combine the 
perturbative calculation with the corresponding power corrections, eventually extending 
the applicability of the result to the peak region. 

As explained in 0, when power accuracy is required, using (^) is cumbersome since 
the number of relevant terms varies as a function of Q and p. To evaluate In J(i/, Q"^) 
we therefore perform explicitly the Borel integral (0). The prescription we use to define 
the integral avoiding the renormalon singularities as well as the technique used for the 
analytic integration were explained in detail in 

To complete the perturbative calculation we need to match the resummed result with 
the exact NLO calculation, which includes non-logarithmic contributions. We do this 
using the so-called logi? matching scheme 0. This amounts to replacing the first two 
terms in the expansion of In R{Q'^,y) (where y = t^pu) with the exact terms in the 
following way: 

\nR{Q\y)l^ = R,iy)a{Q')+[R2{y)-lRiiyr)a\Q') 

2 . NLO 



+ In R{Q\y) ^^^-\nR{Q\y) 



DGE 



Here a = aj^/n, Ri{y) and R2(y) denote the exact LO and NLO coefficients respectively, 

I NLO 

Idge 



and \n R{Q'^ , y)\^^ is the expansion of the DGE resummed result to (9(a^). For Ri{y) 
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we use the analytic results and for R2{y) we use formulae that parametrize the coeffi- 



cients extracted numerically from EVENT2 |40]. Whereas the resummed DGE result is 



independent of the renormalization scal^, the matched result has some scale dependence. 
Ignoring the renormalon contributions which are not logarithmically enhanced is justified 
a posteriori: in P] the residual scale dependence was found to be very small for the thrust 
distribution (< 1% for t < 0.25). We expect a similar dependence for the heavy-jet mass. 

2.3 Limitations of the approximation used 

The perturbative DGE result presented in the previous sections was the basis for the 
phenomenological analysis of the thrust distribution in using (almost) the entire range 
in t. However, as will soon become clear, the range of applicability of the result for the 
heavy-jet mass is more restricted. To understand why, and eventually in what conditions 
a quantitative study of the heavy-jet mass can be performed, we will now return to some 
of the approximations made in the derivation and discuss their validity. 

As any resummation, the DGE result captures specific features of the distribution, 
which show up as large corrections, neglecting other, sub-dominant features. To judge the 
quality of the approximation, we must know what we neglect. The first crucial comparison 
is, of course, with fixed-order results. Although the known fixed-order results (here LO 
and NLO) are taken into account by matching, their comparison with the resummed 
result is a valuable tool to identify the range of parameters where the approximation is 
valid and what it may miss at yet unavailable higher orders (NNLO and beyond). 

The comparison of the 0(0^) term in the resummed result, expanded as 

R{Q', y) = l + Ri{y)a{Q^) + R2{y)a\Q') + ..., 
with the exact NLO term calculated by EVENT2 is shown in Fig. [l|. The difference 



between the two increases for y ^ ?/(peak), particularly because of the non- logarithmic 
terms. The approximation for large values of the variable is much better in the case 
of the thrust compared to that of the heavy-jet mass. For example, the ratio between 
dR!^^^ {t) / dt and dR^^^{t)/dt is larger than 0.8 for t < 0.25, whereas for the heavy-jet 
mass the corresponding ratio reaches 0.8 already at = 0.13. The figure also shows that 
for the thrust distribution a modification of the logarithm to L = ln(l/i: — 1), as suggested 
in gives a better approximation for large t, whereas for intermediate t (~ 0.1) the 
approximation with a modified log is worse. In the case of the heavy-jet mass the effect of 
modifying the log is much smaller than the difference with the exact result. We conclude 
that in the latter case modifying the log cannot be considered as a reliable estimate of 
the uncertainty due to the missing non-logarithmic higher-order terms. We will explain 
why this is so, below. 

Having based our exponentiation on the large /3o limit, rather than on a fixed loga- 
rithmic accuracy, it is interesting to know which logs are computed exactly. As before, 
we categorize the logs as they appear in the exponent. In J(z/, Q^), consistently with 0]: 

*It should be stressed that independence of the renormaUzation scale applies within a given renor- 
malization scheme - the scheme in which as(fc^) is defined. Here the scheme is fixed a la 't Hooft, so 
that the renormalization group equation has only two non-vanishing coefficients. This arbitrary choice 
is reflected in terms that are subleading in /3o at the exponent, at NNLL accuracy and beyond. 
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Figure 1: The NLO coefficient for the thrust (upper set of curves) and heavy-jet mass 
(lower set of curves) distributions obtained by expanding the DGE resuh (dashed and 
dotted hne for the unmodified and modified log, respectively) compared to the exact 
result (solid line). 



leading logs are 0(L"+^q;"), NLL are 0(L'^a"), etc. In addition we must distinguish, at 
the level of the exponent, between terms that are leading in /5o which are directly cal- 
culated (section 2.1) and terms that are subleading in /3o, e.g. 0{CpC a^^^ o?^^ which are 
introduced through subleading terms in the splitting function or through the (two-loop) 
evolution of the couphng. 

Terms that are leading in /?o at the exponent would have been exact if not for non- 



inclusive contributions Such contributions ffist show up due to the emission of a 

soft gluon at large angle, that splits into a quark and an antiquark which end up in op- 
posite hemispheres with respect to the thrust axis. The contribution of such an emission 
to the single jet mass 

In J(z/,Q2) is (^(Cp/^oLa^) i.e. NNLL. However, both the thrust 
and the heavy-jet mass distributions receive only a 0{Cp(5qC?^ correction, i.e. N^^^^LL. 
The reason is that in the latter case, contrary to the former[], the observables and their 
inclusive approximations are of the same order - the order of the soft gluon momentum. 
Consequently, the difference between the integrated cross section and its inclusive ap- 
proximation has support only in a small part of the phase space. The OiCpfi^Lo?^ term 
obtained from expanding the inclusive DGE result indeed coincides with the numerical 
value extracted from EVENT2. This was checked in the case of the thrust in (see the 
appendix). We have now made a similar check for the heavy-jet mass. 



^'We are grateful to Gavin Salam for explaining this point. 
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Figure 2: A contour plot of the doubly differential light and heavy-jet mass cross section 
based on 10 million events generated with Pythia at Q = Mi on the parton level, using 
only light quarks (uds). The contours spacing is based on a logarithmic scale: the lowest 
contour corresponds to two events and the others are spaced by a relative factor of 5. 
This illustrates the available phase space and, in particular, the constraints on pi for 
Ph > 1/6. 

Terms that are subleading in /Sq at the exponent are computed exactly to NLL accu- 
racy for both the thrust (^) and the heavy-jet mass ( pSD distributions, thanks to the use 
of the two-loop renormalization group equation and the gluon bremsstrahlung effective 
charge. The latter includes the full non-Abelian contribution to the singular part of the 
splitting function |^, ^j. As in the Abelian case discussed above, the inclusive approx- 
imation has a price. The effect of a large-angle soft gluon that splits into two gluons 
which end up in opposite hemispheres is larger than in the Abelian case, because the sec- 
ond branching contains an additional soft singularity. Indeed, Dasgupta and Salam have 
shown that the single-jet mass In J(z/, Q^) receives a 0{CfCaL'^C(1) contribution, i.e. 
NLL, from such a configuration. As in the Abelian case, the corresponding contributions 
to the thrust and the heavy jet mass are suppressed, and appear at NNLL accuracy or 
beyond. 

This discussion teaches us that the logarithmically enhanced terms are well under con- 
trol for both the thrust and the heavy-jet mass distributions. The key to understanding 
the difference between these two reflected in Fig. |l| is the kinematic constraints 

on the jet masses, which manifest themselves primarily through non-logarithmic correc- 
tions. 

When a hard gluon is emitted at a moderate angle (the heavy-jet mass is not very 
small) emission of further soft gluons is not controlled by the same matrix element. Also 
the phase space is drastically modified. This is not taken into account in our approx- 
imate calculation, where gluons are emitted off a quark-antiquark dipole. Moreover, 
we explicitly assumed in the derivation of (|28| ) that the masses of the two hemispheres 
are independent. This assumption is consistent with the approximation of independent 
emission we employ in the two-jet region. However, it is clearly violated at larger pn- 

Contrary to the approximation we use, in a realistic phase space, the two hemisphere 
masses are correlated. The easiest way to see this is to consider a slightly different 
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definition of tlie jet masses, wliere tlie liemisplieres are defined by minimizing the sum of 
tlie two jet masses ptot = Pi + P2- In tliis case it was sliown |^ tliat tlie largest possible 
value of Ptot is 1/3. This means that pi cannot exceed min {p//, 1/3 — pi/}, which is 
always smaller than 1/6. This constraint becomes quite stringent for large values of pn- 
It manifests itself through a Sudakov form factor, which resums logarithms of the form 
ln(l/3 — ph)- These large perturbative corrections have nothing to do with the Inpjj 
terms that are resummed in the two-jet region. 

Turning to the jet mass definition based on the thrust axis, the upper limit changes 
to pi + P2 < 1/2, but the value of ph for which pi reaches its maximum is still 1/6 
(this value does not depend on the definition of the thrust axis, since at this point the 
distribution of particles is completely symmetric). The available phase space for pi as 
a function of ph is illustrated in Fig. ^ The figure shows the doubly differential cross 
section based on 10 million events generated with Pythia |^ for Q = Mz at the parton 
level. 

Ignoring the kinematic constraints in the calculation, we definitely overestimate the 
cross section for p^ > 1/6. Since the cross section for ph > 1/6 is not negligible (it is 
of the order of ~ 5% a.t Q = Mz), there is also an impact on the cross section below 
1/6. This explains the significant difference (Fig. |I]) between the exact result and our 
approximation, as pn increases. 

The kinematic constraint on the light-jet mass for p// > 1/6 is much more stringent 
than the general condition that the cross section vanishes at pn = 1/2. Whereas the 
former condition dictates the behaviour of the exact NLO coefficient in Fig. |l], it is the 
latter that is imposed by modifying the log from \nl/pH to ln{l/pH — 1). This explains 
why such a modification does not represent the actual effect of non-logarithmic terms at 
large pH- 

The reason why the effect on the thrust distribution is smaller also becomes clear: 
fixing t at some moderate value does not put any stringent constraint on the light-jet 
mass, nor on any other shape variable. When t reaches 1/3 such constraints appear. 
However, since the cross section for t > 1/3 is small (less than 1% at Q = Mz), the effect 
of ignoring the kinematic constraints on the cross section below 1/3 is insignificant. 

It should also be emphasized that the correlation between the hemisphere masses due 
to the restriction on their sum is negative. This should be contrasted with the dynamical 
correlation discussed in [l7, 18 , which is due to radiation from one hemisphere into the 
other ( "non- inclusive" corrections [^). This correlation is positive ||, |l^, |T8|. Recall 
that the non-inclusive correction at the NLO level was shown p[ (see table 2 there, and 
the explanation that follows) to be small in the case of the average thrust, but it increases 
for higher moments (t^), (t^), etc. Based on the definitions (|I|) and (|), both t and ph 
are smaller in the full calculation than in the inclusive one, in configurations where a 
gluon decays into particles that end up in opposite hemispheres. Therefore, in the full, 
non-inclusive calculation the thrust distribution is somewhat flatter than in the inclusive 
one The same holds for the pn distribution. Thus, ignoring the radiation from one 
hemisphere into the other results in underestimating both the thrust and the heavy-jet 
mass distributions at large t and pn, respectively. On the other hand. Fig. |l| shows that 
ignoring the kinematic constraints as well as the radiation from one hemisphere into the 
other, results in overestimating the cross section at large pu, whereas the effect on the 
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thrust is smaller. We conclude that the dominant feature of the pn distribution missed by 
the current approximation is due to the kinematic constraints (non-logarithmic terms). 
Finally, the main conclusion from this analysis is that, for the heavy-jet mass distri- 



bution, an approximation based only on the two-jet limit (pSj) breaks down already at 
pH = 1/6 and, in general, is not as accurate as the analogous calculation of the thrust 
distribution (0). Of course, the approximation is improved by matching the resummed 
expression with fixed order calculations, but in the absence of a NNLO result, large un- 
certainties are unavoidable. Since the cross section is peaked at rather small pn-, it is 
worth while to attempt a quantitative study of the distribution in the two-jet region, 
using an upper cut on the analysed data at p// = 1/6 or below. In performing such an 
analysis one must keep in mind possible large corrections beyond the NLO. 

2.4 Non-perturbative corrections 

Generally, the effect of hadronization on an event-shape observable can be taken into 
account by power corrections 0^113] • The presence of power corrections can be deduced 



from perturbation theory: the renormalon ambiguity which appears upon resumming the 
perturbative expansion must be compensated on the non-perturbative level. This allows 
us to detect, by perturbative tools, certain non-perturbative corrections. Assuming that 
these corrections dominate, we can parametrize the hadronization corrections based on 
the form of the renormalon ambiguity, using a small set of parameters. 

Since the DGE formulae for the thrust (^) and the heavy-jet mass (^) distributions 
are both expressed in terms of the single-jet mass distribution In J(z/, Q^), it is natural 
to parametrize the power corrections through this quantity. This immediately implies a 
relation between the corrections to the two distributions, a relation that can be confronted 
with data. In parametrizing the non-perturbative corrections through In J(z/, Q^), we will 



effectively regard (|2^) and (^) as non-perturbative relations. This means in particular 
that we will ignore correlations between the hemispheres also on the non-perturbative 
level. This should be contrasted with |T^, where such correlations were found to be 



essential for the consistency between the heavy-jet mass distribution and that of the 
thrust. 

As discussed in the first thing we learn from the perturbative formula for DGE (^) 
is that ambiguities, and therefore power corrections, appear in the exponent. Next, the 
renormalon singularities of In J(z/, Q^) in ( ^5]) (taking into account the factor svo.'ku/ttu 
of Eq. (p^)) appear at half integers u = |, |, | . . ., and at u = 1,2. The singulari- 
ties at M = 1, 2 are due to collinear emission. They indicate power corrections of the 
form ~ l/{Q'^p) and ~ l/{Q^p^)^ which can be neglected in the region pQ / > A. The 
half-integer singularities are associated with the large-angle soft emission. They indi- 
cate P power corrections of the form ~ l/{QpY, where n is an odd integer. These power 
corrections are dominant as p approaches A/Q, and they can be resummed into a non- 
perturbative shape function of a single argument, as originally suggested by Korchemsky 
and Sterman [|13|-||18 . 



Possible modifications of the Borel singularities beyond the large /3o limit were also 
examined in 0. It was shown that, through the modification of the Borel poles into cuts, 
both large-angle, soft gluon corrections of even powers and collinear corrections, whose 
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magnitude is controlled by the same parameters, are present. The phenomenological 
analysis below will follow and include all these power corrections: the large-angle, soft 
gluon corrections of any power will be parametrized by a single argument shape function, 
as done in |]16|, |1^ , and the leading power corrections of coUinear origin will be introduced 
even though the analysis of showed that the effect of the latter is very small. 
The integrated cross section of the single-jet mass is then given by. 



dp 
c 27rii/ 



exp 



z/p + In jP^(z/, Q'') + In /'^(z/, Q 



(31) 



where the sum of the perturbative (PT) and non-perturbative (NP) contributions to the 
exponent is a physical quantity and thus independent of the renormalon regularization 
prescription. Our prescription for the Borel integration, which defines the perturbative 
part In J^"'"(z/, Q^) of Eq. (|23), was described in H]. Note that it is different from the 



cutoff definition used in [0, [T^. 

Thanks to the factorized form in z/-space, Eq. ( pi] ) can be written as a convolution: 



dp R'/{Q'', p - p) /NP(Q^ p) ^ R^iQ'', p) ® /NP(Q^ p) 



where the non-perturbative correction is 

dv 



C 2TTi 



exp 



z/p + ln J^^(z/, Q 



(32) 



(33) 



Given that /np(<5^, p) incorporates the non-perturbative corrections to the single-jet mass 
distribution, and given the relations (^6]) and (pT]), i.e. 



RtiQ\t) = R,{Q\t)(g)R,{Q^t) 
Rp„iQ',PH) = {RpiQ',t)}\ 
the corrections to the thrust and heavy-jet mass distributions can then be expressed as: 

Rt{Q\t) = Rr{Q\t) ® /NP(Q^^) ® fMQ^t) 



RpAQ",ph) = \ RZ^Q'^ph) 



1/2 



/Np(<5^P^^) 



(34) 



Next, to parametrize /np(Q^,p), we use what we learned from the renormalon struc- 
ture of the exponent. Ignoring, for the time being, corrections of coUinear origin, which 
appear beyond the large /3o limit, we start, as in [|^, by writing 



lnJ^P(z/,Q^ 



In J^P(z/A/Q) 



n=l 



(35) 



where we exhibited the dependence of the non-perturbative function J^^{vK/Q) on a 
single argument. According to the renormalon ambiguity pattern of (125|), we expect 
corrections of odd powers but not of even powers. Thus the relevant parameters are A„ 
with odd n. Following |16|, ^ the set of parameters A„ can be traded for a single- argument 
shape function. 



J^P(z.A/Q) = dCfiC) exp {-(uA/Q 



(36) 
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This implies that A„ are the central moments of the shape function, namely, 

POO 

Ai = / dCCfiO 



A, 



dC (C-Aif /(C), 



etc. For small u (large p) it is sufficient to keep the leading term in In J^^iyK/Q), which 
gives 

/(C) = 5(C-Ai). (37) 
Thus, for p ^ A/Q the non-perturbative effects can be approximated by a shift of the 



perturbative distribution, in agreement with [[14| , |23| . Based on the expectation that 



the A2 correction is suppressed, the shift approximation may hold in a rather wide range. 
However, in the distribution peak region higher powers in (|35|) become relevant and a 
shape function is required. For this function we will be using the form suggested in M: 



-biC-biC 



(38) 



which is flexible enough for the second central moment A2 to vanish, as predicted by 
the large Pq ambiguity pattern, but we will usually not impose this condition. Here uq 
is flxed by normalization, J^d(f{() = 1, and g, fci, ^2, ^i, &2 are free parameters to be 
determined by the data. 

Finally, the full non-perturbative function (|33D , which we convolute with the pertur- 
bative distribution (ETj, is 



/np (Q^p) 



A V A 



S{p) - A2 



A^ • 

Wp 



(39) 



where the first brackets contain the main, large-angle, soft gluon correction ( pHf ) and the 
second brackets incorporate the leading coUinear correction (see [^]). 

Another aspect that should be taken into account is that the perturbative calculation 
of the distributions is for massless partons, whereas the detected hadrons are massive. 
The effect of finite hadron masses was studied in detail for average values of event-shape 
variables by Salam and Wicke [^. They define two types of corrections: universal, 
which arise from the reshuffling of momenta when forming massive hadrons from massless 
partons, and non-universal, which depends on the way the masses are treated in the 
deflnition of the variables. Based on kinematic arguments the mass effects were found 
to be of order A/Q, whereas a more detailed QCD analysis, using local parton - hadron 
duality, indicated that both corrections scale as (lnQ/A)^A/Q, where A = Ca/Po ~ 1-6. 
The origin of this logarithmic dependence is the fact that the mass effects associated 
with transforming kinetic energy into hadron masses grow with particle multiplicity, 
which increases logarithmically with the energy. 

The comparison between the hadronization corrections to the thrust and those to 
the heavy-jet mass is particularly problematic because the thrust is deflned using three- 
momenta, whereas the heavy-jet mass using four-momenta. According to this implies 
that the two variables will have different sensitivity to the hadron masses, i.e. different 
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non- universal mass effects. This problem can be by-passed |^|, modifying the definition 



for the variables. In addition to the actual definition, the "hadron level" at which the 
measurement is performed makes a difference: depending on the experimental set-up, 
different particles may effectively be stable. 

The non-perturbative parameters used to describe the hadronization corrections may 
be defined in different ways. In particular, they depend on the specific definition of 
the variables, on the "hadron level" used for the measurement, and on the way non- 
perturbative corrections are separated from the perturbative sum. Moreover, there is 
no unique way to distinguish between the hadronization process (the "standard" power 
corrections) and the formation of massive particles (the "mass effects"). Therefore, con- 
trary to as, where the dependence on the particular procedure must be interpreted as 
some "theoretical uncertainty", the numerical values of the non-perturbative parameters 
are meaningful only within a given procedure. From this point of view it is satisfactory 
to have one procedure by which non-perturbative parameters to different event-shape 
observables can be compared. 



In this paper we will consider the following procedures |2J] to define the variables, 
given that the detected particles are massive. We shall refer to these procedures as hadron 
mass schemes (HMS): 

• Massive scheme (M): Use the massless definitions as given in Eqs. (|I])-(^, ig- 
noring the fact that the detected hadrons are massive. This is the scheme that is 
normally used for the published data. 

• P scheme: Use the measured three-momenta Pi as they are and assign Ei = \pi\ 
for each particle. This way three-momentum is conserved but not energy. 

• E scheme: Use the measured energies Ei as they are but rescale the three-momenta 
with a factor Ei/\pi\. This way energy is conserved but not three-momentum. 

• Decay scheme (D): Let all particles (including protons) decay to massless parti- 
cles. In the Monte Carlo this is done by activating the decays of all particles that 
are normally considered stable (/x^, vr^, K"^, K^) . In addition, we let the protons 
(and neutrons for simplicity) decay semileptonically according to p mev, where 
m is massless. This yields a final state consisting entirely of massless particles (the 
electron mass can safely be neglected). This way both energy and momentum are 
conserved, but on the other hand the additional decays are not due to the strong 
interaction. 

For the thrust and heavy-jet mass the non-universal mass effects that are present in the 
conventionally used M scheme are diminished (or become common to the two variables) 



in the decay, P and E schemes [^. The decay scheme, in which energy and momen- 
tum conservation hold and all "measured" particles are strictly massless, seems the most 
attractive for comparison with perturbative calculations. This will therefore be our pre- 
ferred scheme. This choice is natural also according to the picture that emerges from the 



Ariadne |^ Monte Carlo based study in (see Fig. 8 there): the logarithmic depen- 
dence of the power correction coefficient is the smallest in the decay level with respect to 
other hadronic levels. 
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To transform the data from one HMS to another we use the Pythia Monte Carlo 
model (version 6.158). We have verified that the model describes well the thrust and the 
heavy-jet mass distributions in the entire energy range and in a wide range of the shape 
variables. Comparison with thrust data gives a of 154 for 225 points for 14 < Q < 
189 GeV and 0.01 Mz/Q <t < 0.30. Comparison with the heavy-jet mass data gives a 

of 129 for 209 points in the same ranges. 

Using Pythia we calculate transfer matrices T such that the cross section in bin i in 
scheme A is given by af = J2j T^aj., where aj is the cross section for bin j in the original 
HMS. The matrices T are calculated for each HMS, energy and experiment separately. 
For the data at Mz we use 10^ events to calculate the matrices, whereas for the other 
energies we use 10^ events. This way the statistical error in the calculation of T can 
be neglected with respect to the statistical error in the data. It should be emphasized 
that this procedure is very similar to what is normally used to make corrections to the 
data (e.g. those due to imperfect detector efficiency). To estimate the systematic error 
the calculation of T one could be repeated using a different Monte Carlo model such as 



Herwig 1 44]. However, to do it properly one should use the same Monte Carlo for the 



complete analysis chain used for extracting the data. 



3 Data analysis 

In this section we compare the calculations of the thrust and the heavy-jet mass distri- 
butions to the world data and perform fits to extract and the parameters that control 
the hadronization corrections. A similar analysis for the thrust data was performed in 
Here, our main motivation is the comparison between the thrust and the heavy-jet mass. 
To make a fully quantitative comparison we must repeat the thrust analysis taking into 
account several experimental aspects that were not fully considered before. This includes 
the following issues: 

a) Using data where the measurement is corrected to be a measurement of 
all particles. Most experiments base their analysis on charged tracks and then cor- 
rect to include all particles, using a Monte Carlo model such as Pythia. However, 
some experiments have chosen to present data for charged particles only. While 
this procedure gives better control of the systematic error, the measured distribu- 
tion does not correspond to the distribution calculated perturbatively. Thus, we 
want to use data that have been corrected to represent all particles. Out of the 
data sets we used in our earlier analysis the important ones, which are for charged 
particles only, are the ALEPH and DELPHI data sets at Mi- The difference in the 
"ideal measurement" used by different experiments to correct the data explains the 
discrepancy we found between the different data sets at this energy. In order to 
be able to use data from all experiments, especially at Mz, we apply the method 



described in section |2.4| to transform the charged particle data to be for all particles. 



We have not tried to estimate the systematic errors in this procedure. 

b) Using a proper HMS, common to the two variables. The published data for 
the thrust and the heavy-jet mass is based on the definitions (0) and (||), respec- 
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tively. While for the former the definition coincides with the P scheme, the latter 



is strictly a "massive scheme". Following and our discussion in section 
comparison of the hadronization corrections of the two distributions requires to 
use a common, non-massive scheme. Our main analysis will be made in the decay 
scheme. However, we will also compare the different HMS and verify that as ex- 
tracted from the fits does not change (recall that a strong dependence was found 
in 12 



c) Estimating the effect of heavy primary quarks. The published data are 
based on all events, including those where heavy quarks are produced, while the 
calculation always refers to 5 light flavours. The contribution of the quark mass to 
the variable, as well as the difference in the soft-gluon radiation from heavy quarks 
compared to that from light quarks can be a source of error. In order to estimate 
this error, we use Pythia to calculate the distribution based on all (udscb) or only 
light (uds) primary quarks and then use the ratio Ri of the two distributions to 
transform the data bin by bin according to. 

Similarly to the HMS transformations, the calculations of the Ri are done for each 
HMS, energy and experiment separately, taking into account the different binnings 
of the data. 

Some comments are due concerning the data we use. As in 0, H, the dat a are in the 
range 14 to 189 GeV. When several data sets are available from the same experiment and 
energy, we use the latest one. We have removed some data sets with large errors, and 
restricted ourselves to published data, where the details on the Monte Carlo corrections 
(in particular, on whether the data correspond to all particles or to charged particles 
only) are available. The data used for the fits are summarized in Table ^ below. When 
available separately, the systematic errors have been added in quadrature to the statistical 
ones. When fitting the data we have calculated the integrated cross section for each bin 
included in the fit. However, in the plots showing the data we have followed the customary 
practice of plotting the data points at the bin centres, even though this may be somewhat 
misleading for large bin widths. 

Given the theoretical uncertainties for the heavy-jet mass distribution we start with a 
separate analysis of the thrust data, then continue with a separate analysis of the heavy- 
jet mass data, and finally analyse the heavy-jet mass data using the results of the thrust 
analysis as input. This makes it possible to test the DGE method and the assumptions 
we make on the hadronization corrections. 



3.1 Thrust distribution 

We present here the analysis of the thrust distribution. Similarly to we begin by per- 
forming fits based on a shift of the DGE perturbative distribution in a limited range 



and then turn to shape-function-based fits |]l6i \[% in a wider range. As explained above, 
we will also study the HMS dependence and the effect of heavy quarks as well as the 
impact of using data for all particles instead of charged particles only. 
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Figure 3: Thrust data at Mz in the peak region together with shift-based fits using the 
range 0.05 Mz/Q < t < 0.30 (dashed hne) in different HMS. The perturbative calculation 
(DGE) is also shown (dotted line). 



3.1.1 Fits with a shift 

The analysis is based on fitting as and a shift \iK/Q of the perturbative distribution 
according to (pT]). This applies to t values far enough above the peak region, K/Q. 

We start by investigating the dependence on the HMS and the effect of heavy primary 
quarks. Figure |] shows the fit results together with the data at Mi in the different HMS. 
We see that both the position and the height of the peak vary. Note also the agreement 
between the different data sets in the peak region in a given HMS (cf. Fig. 9 in 0) when 
the ALEPH and DELPHI data have been corrected from charged particles only to all 
particles. 

The results of fitting and Ai to the data are summarized^] in Table |l]. The global 

"''When comparing with |9|, the reader should keep in mind that here the non-perturbative correc- 
tions correspond to the single-jet mass lnJ(z/, Q^) in ( [2^ ) rather than to the exponent of the thrust, 
2 In J(z^, Q^). Thus, the moments Ai, A2, etc., are smaller by a factor of 2. 
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Table 1: Fits of as and Ai to the thrust data with a lower cut tmm = 0.05 Mz/Q and an 
upper cut tmax = 0.30 in different HMS. 



HMS 


aMs(^z) 


AiA [GeV] 


xVdof 


Points 


M/P (udscb) 


0.1072 ± 0.0004 


0.35 ± 0.02 


0.58 


159 


M/P (uds) 


0.1074 ± 0.0005 


0.31 ± 0.02 


0.58 


159 


E (udscb) 


0.1080 ± 0.0004 


0.42 ± 0.02 


0.51 


159 


E (uds) 


0.1082 ± 0.0004 


0.38 ± 0.02 


0.51 


159 


D (udscb) 


0.1086 ± 0.0004 


0.63 ± 0.02 


0.43 


159 


D (uds) 


0.1088 ± 0.0004 


0.59 ± 0.02 


0.45 


159 



Table 2: Fits of as and Ai to the thrust data in the decay scheme (udscb) as a function 
of the lower cut tmin = i^mm(^z) Mz/Q with the upper cut fixed as tmax = 0.30. 



tmin(Mz) 




ays{Mz) 


AiA [GeV] 


xVdof 


Points 


0.03 


0. 


,1075 ± 0.0002 


0.69 ± 0.01 


0.55 


190 


0.04 


0. 


,1075 ± 0.0003 


0.69 ± 0.01 


0.54 


173 


0.05 


0. 


,1086 ± 0.0004 


0.63 ± 0.02 


0.43 


159 


0.06 


0. 


,1091 ± 0.0005 


0.60 ± 0.03 


0.42 


147 


0.07 


0. 


,1091 ± 0.0007 


0.60 ± 0.04 


0.45 


135 


0.08 


0. 


,1090 ± 0.0007 


0.61 ± 0.04 


0.43 


124 


0.09 


0. 


,1086 ± 0.0009 


0.64 ± 0.06 


0.45 


111 


0.10 


0. 


,1084 ± 0.0011 


0.65 ± 0.07 


0.48 


104 



impact of ignoring the difference between heavy and light primary quarks on the fit is 
very small and it will be neglected in most of the following analysis. This does not imply, 
of course, that the effect is small at low energies. It is quite clear that it is significant at 
14 GeV, but the weight of these data in the fit is small. The table also shows that different 
HMS yield very similar values of (the differences are less than 1.5%), whereas Ai 
changes by a factor of ~ 2. The stability of as is very important: a physical parameter 
must not depend on arbitrary choices such as the particular definition of the variable or 
the "hadron level" . Note that stability was not reached in , when extracting as from 
mean values of event-shape variables. Analysing the distribution is advantageous, since 
the dependence on the variable, in addition to the energy, constrains ag further. We have 
explicitly checked, performing a NLL fitj^ (with log R matching) to the distribution that 
the results show the same stability as in the DGE case. From the table it is also clear 
that the decay scheme gives the best fits. We will mainly use this HMS in the following. 

Next, we consider the dependence on the lower cut of the fit, tmin. The results are 
summarized in Table |^. Note that for t^^-^ > 0.05 Mz/Q the results are stable: all fits 
for larger tmin agree within errors with the result for tmin = 0.05 Mi/Q- This means that 
there is essentially no dependence on the lower cut (it is less than 0.5%), provided that 

§The result in the decay scheme for all primary quarks is a-^{Mz) = 0.1217 ± 0.0006 and AiA = 
0.69 ± 0.02 GeV with xVdof = 0.42. 
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the cut is far enough above the peak region. 

We have checked that varying the upper cut tmax between 0.15 and 0.35, while keeping 
the lower cut tmin = 0.05 Mz/Q fixed, the results for aMs{Mz) and Ai agree, within errors, 
with the results obtained for tmax = 0.30. Thus, the dependence on the upper cut is also 
negligible (less than 0.5%). 

To demonstrate the importance of using data for all particles rather than for charged 
particles only, we made a separate fit to the ALEPH and DELPHI data at Mz in the P 
scheme. The result changes from aMs(Mz) = 0.1079 ± 0.0005 and AiA = 0.33 ±0.03 GeV 
with xVdof = 0.45 for all particles, to aMs(Mz) = 0.1091 ±0.0005 and AiA = 0.27 ±0.03 
GeV with x^/dof = 0.62 for charged particles. Thus, the difference in the fit results is 
small, but not negligible. However, there is a significant difference in the quality of the 
fit: in addition to the lower x^/dof, the sensitivity to varying the cuts on the fitting range 
is much lower. This concerns in particular the lower cut - compare Table Q above with 
Table 2 in 0. 

Finally, knowing that beyond the large-/?o limit power corrections are modified by 
logarithmic Q dependence |P] , it is interesting to look for logarithmic Q dependence of Ai 
by performing separate fits at each energy. Additional motivation for such investigation 
is provided by the finding of that the scale dependence of hadron-mass related non- 



perturbative effects is enhanced by a factor of \n'^{Q/A), where A = 1.565. Fixing as to 
the value obtained by fitting all energies, we fitted Ai separately for each energy. The 
results of these fits in the decay and P schemes are shown in Fig. ^ along with a fitted 
line of the form AiA = pi + P2 ln(Q/A). The best-fit values of pi and p2 are 

^ , . , , Pi =0.13 ±0.11 GeV 

D (udscb) : ^ Q ^ Q GeV, ^^^^ 

PMscb)- Pi = -0.16 ±0.12 GeV 

r (udscDj . ^ ^ Q GeV. ^^^^ 

The results show some trend of increase in Ai, but the power of the logarithmic de- 
pendence cannot be extracted from the data. The difference between the decay and P 
schemes appears to be large, but independent of Q. In the following we will ignore the 
logarithmic Q dependence and, instead, concentrate on fits with a shape function and 
the determination of A2. 

To summarize, the best fit values of and Ai when using the decay scheme with all 
primary quarks, are 

/ 1 UN aMs(Mz) = 0.1086 ±0.0004 
D (udscb): ;^X= 0.63 ±0.02 GeV ^^^^ 

with x^/dof = 0.43. These values were found to be stable under variations of the lower 
and upper cuts of the fitting range and to have very small uncertainty due to heavy 
primary quarks. The uncertainty in the determination of ag based on varying the HMS 
was found to be ~ 2%. The shift Ai can be a meaningful parameter only when the HMS 
has been fixed, in which case it is determined up to 5 - 10% depending on the fitting 
range and the effect of heavy primary quarks. In the different schemes we considered, Ai 
changes by up to a factor ~ 2, giving an indication on the significance of hadron mass 
effects. 
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Figure 4: The Q-dependence of the shift Ai for a fixed in the decay and P schemes. 
The hnes are fits of the form AiA = pi + P2 ln(Q/A). 



3.1.2 Fits with a shape function 

Fitting with a shape function has the advantage that the peak region can be included in 
the analysis and the structure of subleading power corrections can be detected. While the 
motivation is to extend the fitting range as much as possible, the region t ~ A/Q, where 
all powers are equally important, must be avoided. Similarly to ||^, these considerations 
lead us to choose t^i^ = 0.01 Mz/Q — lA/Q. The following fits are made to data in the 
decay scheme (with all primary quarks). 

For the functional form of the shape function (pSD, we use two different possibilities 
for the fall-off at large (: exponential and Gaussian. The results for an exponential fall-off 
are: 

a^(Mz) = 0.1090 ± 0.0004 

AiA = 0.605 ± 0.013 GeV AaA^ = 0.002 ± 0.024 GeV^ ,,,, 
6, = 0.61 ±0.10 62 = (fixed) ^^^^ 

q = 1.55 ±0.11 

with x^/dof = 0.53 for 225 points, and those for a Gaussian fall-off: 

aMs(^z) = 0.1085 ± 0.0004 
^ , J AiA = 0.625 ±0.011 GeV AaA^ = -0.049 ± 0.018 GeV^ 

D (udscb) : ^ Q ^^^^^^ ^ Q_Q28 ± 0.004 ^^^^ 

q = 0.77 ±0.08 

with x^/dof = 0.58 for 225 points. Both results are in good agreement with the shift- 
based fits for tmin > 0.05 Mz/Q (^3]). As was the case in our earlier analysis, the best fit 
turns out to have an exponential fall-off. The Gaussian fall-off gives a slightly worse x^- 
In addition, the softer, exponential fall-off is closer to the asymptotic behaviour found 



in |18| in a calculation of the shape function using a perturbative model. 
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Table 3: Fits of aMs(Mz) and shape function with exponential fall-off to the thrust data 
in the decay scheme (udscb) as a function of the lower cut tmin = ^min(^z)A^z/Q with 
the upper cut t^ax = 0.30 fixed. 



tmin(Mz) 


aMs(Mz) 


AiA [GeV] 


A2A' [GeV^] 


xVdof 


Points 


0.01 


0.1090± 0.0004 


0.605 ± 0.013 


0.002 ± 0.024 


0.53 


225 


0.02 


0.1088± 0.0004 


0.619 ± 0.015 


-0.004 ± 0.023 


0.48 


209 


0.03 


0.1089± 0.0004 


0.618 ± 0.016 


-0.001 ± 0.023 


0.47 


190 


0.04 


0.1089± 0.0006 


0.623 ± 0.040 


0.015 ± 0.093 


0.42 


173 



The sensitivity to the functional form gives a good indication of the difficulty in 
determining A2 and higher moments of the shape function from the data in the absence 
of more definite theoretical input. The results of the exponential fall-off fit, which is 
the best fit, are in agreement with the large (3o renormalon ambiguity pattern, which 
suggests that A2 = 0, but this is not true for the Gaussian fall-off. Nevertheless, in 
both cases A2 -C A^ Thus, if AiA is used as a reference scale for the non-perturbative 
parameters, then A2 is always small. One should also keep in mind that the value of A2 
is strongly correlated to the values of and Ai. For example, fixing the values of and 
Ai according to the fit with an exponential fall-off when fitting with a Gaussian fall-off 
the result is A2A^ = —0.020 ± 0.009 GeV^, which agrees, within errors, with the result 
for a exponential fall-off. 

We have also verified that the result of the shape-function-based fit is not sensitive 
to the fitting range. This is true for as as well as for the non-perturbative parameters Ai 
and A2. In particular, the dependence on the lower cut tmin is summarized in Table ^ It 
is interesting to note that the fact that A2 is small can be deduced even if the region to 
the left of the distribution peak is excluded from the fit. 

Figure | shows the data at Mz in the decay scheme together with the best shape- 
function-based fit ( PI) and a shift-based fit. For the latter we use the range 0.06 Mz/Q < 
t < 0.30, which gives the same as as the shape-function fit (see Table H). The figure illus- 
trates clearly the hierarchy between power corrections as a function of t. For large t the 
shape function is well approximated by a shift. In the peak region the non-perturbative 
corrections with higher power become important, but the effects are still small. Finally, 
below the peak all power corrections become equally important. 

For illustration, the best shape-function-based fit (H ) is also shown at all energies, 
together with the data, in Fig. |^. The contribution to the from the individual experi- 
ments for the best shape-function-based fit (0) and the shift-based fit is given in Table ^ 
The fits are good at all energies, except at 14 GeV. This exception is not surprising, since 
the effects of heavy quarks were not taken into account. 

To investigate the effects of heavy primary quarks we have also performed the shape- 
function-based fit to the data which was transformed according to (^OD , so as to remove 
heavy-quark events from the data sample. As before, this transformation has been made 
using the Pythia Monte Carlo model. The results using a shape function with exponen- 
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Figure 5: Thrust data at Mz in the decay scheme together with the best shape-function- 
based fit (sohd hne) and a shift-based fit using the range 0.06 Mz/Q <t < 0.30 (dashed 
hne). The DGE perturbative calculation is also shown (dotted line). 



tial fall-off, are 



ayd.Mz) = 0. 1092 ± 0.0004 

AiA = 0.568 ± 0.013 GeV AaA^ = 0.005 ± 0.023 GeV^ 



6i=0.64±0.10' 62 = ■ ■ 

q = 1.51 ±0.12 

with x^/dof = 0.63 for 225 points. This shows (cf. Table |lD that the effect of heavy 
quarks is slightly larger in the peak region with respect to higher t, but it is still small. 

Considering the contribution to the from the individual experiments, that of 
14 GeV becomes only slightly smaller (14.1 12.2) whereas all the others increase. 
The Monte Carlo results suggest that the effect of heavy quarks is small, except at 14 
GeV, in agreement with [31]. However, the change in raises doubts concerning the 
reliability of using the Monte Carlo together with Eq. ( pO]) to quantify this effect. It is 
clear that there is place for improving the analysis in this respect. 

Finally, we give the results of a fit with a shape function with exponential fall-off to 
the data in the P scheme: 

aMs(^z) = 0.1068 ± 0.0004 
^ , J AiA = 0.388 ±0.013 GeV AsA^ = 0.069 ± 0.022 GeV^ 

P(udscb): ^^^0.47 ±0.05 h^O ^^^^ 

q = 0.73 ±0.08 

with x^/dof = 0.67 for 225 points. This result agrees reasonably well with the shift-based 
fit in Table |I|. The fact that here, in contrast to (^^, the non-perturbative parameter 
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Figure 6: Comparison of the thrust data in the decay scheme with the best fit shape 
function. Both the fit and the data have been multiphed by a factor 20^, with = . . . 9 
for Q = 14 . . . 189 GeV. Data points with errors larger than 100% are not shown. 



A2 is significantly different from zero, A2 ~ A^, is interesting. The differences between 
the non-perturbative parameter^ A^ in different HMS are associated with mass effects. 
It is rather clear that these effects, at least their non-universal parts do not have 
a "perturbative" origin, and therefore it is not surprising that they do not admit the 
renormalon power correction pattern. 

To summarize the fits with a shape function, the best fit in the decay scheme is (0), 
which agrees well with the shift-based fit (|43|) . The second central moment A2 in this 
scheme is consistent with zero, in accordance with the large /3o renormalon pattern. This 
is not the case in other HMS. Similarly to the shift-based fits, the only non- negligible 
source of uncertainty in as is the choice of HMS. 

^It is expected that the finite hadron mass aflFects also the higher power corrections. Note that the 
modification of the distribution when changing the HMS is more comphcated than just a shift - see 
Fig. I 
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Table 4: Summary of tfie experimental data used and the contribution to from the 
individual experiments (in the decay scheme with all primary quarks). For the thrust, 
the data were fitted using a shift (0.06 Mi/Q <t< 0.30) and a shape function (SF) with 
exponential fall-off (0.01 Mi/Q <t< 0.30), whereas for the heavy-jet mass (0.01 Mi/Q < 
Ph < 0.10) the values were obtained with aMs(Mz) = 0.109 and with the shape 
function fixed according to the thrust analysis. 













thrust 


shift 


thrust 


SF 


heavy jet SF 
(a, & SF fixed) 


Experiment Reference Q [GeV] 




Points 




Points 


x' 


Points 


TASSO 




16 


14 


0.00 





14.12 


6 


1.59 


1 


TASSO 


45, 


16 


22 


0.00 





3.80 


7 


0.52 


2 


JADE 




47 




35 


2.16 


5 


8.59 


11 


8.60 


3 


TASSO 


45, 


16 


35 


1.44 


3 


3.94 


8 


0.61 


3 


JADE 




m 


44 


4.24 


6 


4.82 


11 


4.82 


3 


TASSO 


45, 


16 


44 


0.84 


3 


4.22 


8 


0.08 


3 


ALEPH 




48| 


91.2 


3.35 


9 


7.82 


17 


12.10 


10 


DELPHI 




49 




91.2 


9.60 


12 


14.64 


17 


10.28 


8 


L3 




5C| 


91.2 


3.67 


7 


3.90 


9 


6.15 


5 


OPAL 




5T 


1 


91.2 


10.90 


24 


20.42 


29 


66.23 


10 


SLD 




52 




91.2 


3.48 


5 


5.19 


7 


11.69 


2 


ALEPH 




53 




133 


1.58 


6 


1.60 


8 






DELPHI 




54 




133 


1.47 


5 


1.87 


7 


0.62 


3 


OPAL 




55 


1 


133 


2.16 


6 


2.51 


10 


2.09 


6 


DELPHI 




54 


1 


161 


3.47 


6 


4.81 


7 


1.06 


3 


OPAL 




56] 


161 


1.77 


7 


1.95 


10 


1.70 


6 


DELPHI 




54 


1 


172 


3.11 


6 


3.18 


7 


0.18 


3 


OPAL 




57 


1 


172 


2.14 


7 


2.31 


10 


2.74 


7 


DELPHI 




54 




183 


4.22 


14 


4.84 


16 


0.62 


7 


OPAL 




57 


1 


183 


0.58 


8 


0.62 


10 


2.04 


7 


OPAL 




1 


189 


0.94 


8 


1.40 


10 


2.32 


7 


Sum 


61.1 


147 


116.5 


225 


136.0 


99 



3.1.3 Fits with a modified logarithm 

To estimate the effects of the non-logarithmic terms, which are not included in the re- 
summation, we repeat the calculation using a modified logarithm, L = ln(l/t — 1) 0. 
The results of fitting as and a shift Ai using the modified logarithm in different HMS 
is summarized in Table ^. Comparing with Table |l] we see that the central values of as 
are ~ 2% larger and AiA is ~ 0.1 GeV smaller, whereas the HMS dependence is un- 
changed. We also note that the values are significantly worse when using a modified 
logarithm. 

In addition to the HMS dependence we have also investigated the dependence on the 
lower and upper cuts for the fit in the decay scheme. Recall that with the unmodified 
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Table 5: Fits of as and Ai to the thrust data with a lower cut tmin = 0.05 Mz/Q and 
an upper cut tmax = 0.30 in different hadron mass schemes using a modified logarithm, 
L = ln(l/t - ij^ 



HMS 


aMs(Mz) 


AiA [GeV] 


xVdof 


Points 


M/P (udscb) 


0.1096 ± 0.0005 


0.25 ± 0.03 


0.70 


159 


E (udscb) 


0.1104 ± 0.0005 


0.32 ± 0.02 


0.62 


159 


D (udscb) 


0.1107 ± 0.0005 


0.55 ± 0.02 


0.50 


159 



Table 6: Modified logarithm (L = ln(l/t — 1)) fits of and Ai to the thrust data in 
the decay scheme (udscb) as a function of the lower cut tmin = im\n{,Mz)Mz/Q with the 
upper cut fixed as tmax = 0.30. 



tmin(Mz) 


aMs(^z) 


AiA [GeV] 


xVdof 


Points 


0.03 


0.1086 ± 0.0003 


0.65 ± 0.01 


0.73 


190 


0.04 


0.1090 ± 0.0003 


0.64 ± 0.01 


0.72 


173 


0.05 


0.1107 ± 0.0005 


0.55 ± 0.02 


0.50 


159 


0.06 


0.1119 ± 0.0006 


0.48 ± 0.03 


0.44 


147 


0.07 


0.1129 ± 0.0008 


0.42 ± 0.04 


0.45 


135 


0.08 


0.1134 ± 0.0009 


0.39 ± 0.05 


0.44 


124 


0.09 


0.1138 ± 0.0012 


0.36 ± 0.08 


0.46 


111 


0.10 


0.1146 ± 0.0016 


0.31 ± 0.10 


0.48 


104 



logarithm we found no dependence on these cuts outside the peak region. However, 
with the modified logarithm we find a large tmin dependence while keeping the upper 
cut fixed at tmax = 0.30. As shown in Table ^ changing the lower cut from 0.05 Mz/Q 
to O.lOMz/Q, the result for as increases by 4% and AiA decreases by ~ 0.2 GeV. The 
values only stabilize for tmm = 0.07 Mz/Q, giving a-y^g^Mz) = 0.113 ± 0.001 and AiA = 
0.42 ± 0.04. Likewise we find that keeping the lower cut fixed at tmin = 0.07 Mz/Q and 
varying the upper cut from 0.15 to 0.35 the result for as changes by 1.5% and AiA by 
~ 0.05 GeV. 

Even though the fits with a modified logarithm are significantly worse than the fits 
with an unmodified logarithm these results can be used to estimate the significance of 
the non-logarithmic terms. Based on this we conclude that the effects of non-logarithmic 
terms can be as large as 4% in (from aj^{Mz) = 0.109 to 0.113). This is by far 
the dominant contribution to theoretical uncertainty. We emphasize that this is a rough 
estimate and full NNLO calculations are essential to reduce the uncertainty. 

3.2 Heavy jet mass distribution 

As explained in section |2.3| , the approximation we use for the phase space breaks down 
completely at = 1/6. Contrary to the case of the thrust, non-logarithmic corrections 
reflecting additional kinematic constraints may be important at rather small values of 
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Table 7: Fits of and Ai to the heavy-jet mass data with a lower cut of p, 
0.03 Mz/Q and an upper cut of pmax = 0.15 in different HMS. 



HMS 


ttMs(A^z) 


AiA [GeV] 


xVdof 


Points 


P (udscb) 


0.1009 ± 0.0006 


0.58 ± 0.03 


0.42 


87 


P (uds) 


0.1012 ± 0.0008 


0.55 ± 0.04 


0.41 


87 


E (udscb) 


0.1014 ± 0.0007 


0.69 ± 0.04 


0.44 


87 


E (uds) 


0.1018 ± 0.0007 


0.65 ± 0.04 


0.43 


87 


D (udscb) 


0.1019 ± 0.0007 


0.88 ± 0.03 


0.50 


87 


D (uds) 


0.1020 ± 0.0007 


0.85 ± 0.03 


0.50 


87 


M (udscb) 


0.1026 ± 0.0006 


1.03 ± 0.03 


0.56 


87 



Ph- This will be reflected in the independent phenomenological analysis of the heavy-jet 
mass distribution in section 3.2.1. We prefer to rely on the thrust analysis as much as 
possible and compare between the non-perturbative corrections to the two observables 
only in the limited region where the calculation of the heavy-jet mass distribution can 
be trusted. This is the purpose of section 3.2.2. 

3.2.1 Independent fits 

The purpose of this section is to confront the DGE resummed cross section and the 
suggested parameterization of power corrections with the heavy-jet mass distribution 
data and find how severe the limitations described in section |2.3| are in practice. 

We begin by fitting as and a shift of the perturbative distribution (pT]). This should 
be applicable far enough above the peak region. However, since the region where the 
calculation can be trusted is strongly restricted from above, we end up with a rather 
narrow range. The upper cut will be chosen as pmax = 0.15. We will show that the lower 
cut can be as small as pmin = 0.03 Mz/Q. 

The dependence on the HMS and the effect of heavy primary quarks are summarized 
in Table |^. First of all we note that the extracted values of as are much smaller than 
for the thrust case. The reason is that our approximation tends to overestimate the 
higher order perturbative coefficients at large pn'- an example is provided by the NLO 
coefficient shown in Fig. ^ This statement will be further supported by studying the 
dependence on the fitting range (Table It should also be kept in mind that, since the 
distribution is normalized, overestimating the cross section at large pn is compensated by 
underestimating it at low pn- For the NLO coefficient the turnover point is pn ~ 0.05. 

Since the values of as obtained here are so different from the thrust case, it is no 
surprise that the values of the shift Ai differ by much as well. In any case, a meaningful 
comparison of these parameters requires the same value of as to be used. We will come 
back to such a comparison below. 

We note that, similarly to the thrust case, the HMS dependence and the effects of 
primary heavy quarks on the value of as is small. The largest difference is between the P 
and the massive schemes. The differences between the two can be clearly seen in Fig. ^ 
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Figure 7: Heavy jet mass data at Mz in the peak region together with shift-based fits 
using the range 0.03 Mz/Q < pn < 0.15 (dashed hne) in different HMS. The perturbative 
calculation is also shown (dotted line). 



where the fit results at Mz in the different HMS are shown together with the dataU]. 
In the massive scheme the contribution from the hadron masses themselves shifts the 
distribution to the right. In general, as in the thrust case, the HMS transformations 
modify both the position and the height of the peak. Note, however, that the decay and 
massive schemes are quite similar for the heavy-jet mass. The reason is that the only 
difference between the two is the definition of the thrust axis. 

Next, we investigate the variation of the result as a function of the fitting range in 
the decay scheme. From Tables |] and ||, it can be concluded, based on the same criteria 
applied in the thrust analysis, that the maximal range that can be used is 0.03 Mz/Q < 
Ph < 0.15. However, contrary to the thrust case, here we identify a clear trend: as the 
fitting range shifts towards larger values of pn, the value of as decreases. Thus, there is 

"The figure also illustrates the agreement between the different data sets in the peak region (in a 
given HMS) once they are all corrected to be for all particles. 
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Table 8: Fits of as and Ai to the heavy-jet mass data in the decay scheme (udscb) as a 
function of the lower cut pmin = pnim(^z)^z/<55 with the upper cut fixed as pmax = 0.15. 



Pmm(Mz) 




AiA [GeV] 


xVdof 


Points 


0.02 


0.1009 ± 0.0005 


0.94 ± 0.02 


0.62 


106 


0.03 


0.1019 ± 0.0007 


0.88 ± 0.03 


0.50 


87 


0.04 


0.1023 ± 0.0010 


0.86 ± 0.05 


0.44 


72 


0.05 


0.1012 ± 0.0014 


0.95 ± 0.09 


0.38 


60 


0.06 


0.1000 ± 0.0019 


1.05 ± 0.14 


0.34 


48 


0.07 


0.0990 ± 0.0031 


1.12 ± 0.24 


0.36 


43 



Table 9: Fits of as and Ai to the heavy-jet mass data in the decay scheme (udscb) as a 
function of the upper cut pmax, with the lower cut fixed as pmin = 0.03 Mz/Q. 



Pmax 


aMs(^z) 


AiA [GeV] 


xVdof 


Points 


0.10 


0.1028 ± 0.0009 


0.86 ± 0.04 


0.31 


58 


0.11 


0.1027 ± 0.0009 


0.87 ± 0.04 


0.30 


60 


0.12 


0.1024 ± 0.0008 


0.87 ± 0.03 


0.38 


69 


0.13 


0.1023 ± 0.0008 


0.87 ± 0.03 


0.40 


76 


0.14 


0.1021 ± 0.0007 


0.88 ± 0.03 


0.43 


84 


0.15 


0.1019 ± 0.0007 


0.88 ± 0.03 


0.50 


87 


0.16 


0.1014 ± 0.0006 


0.90 ± 0.03 


0.57 


95 


0.18 


0.1006 ± 0.0006 


0.92 ± 0.03 


0.93 


103 


0.20 


0.0995 ± 0.0005 


0.95 ± 0.03 


1.39 


114 



a systematic discrepancy between the shape of the theoretical distribution and the data. 
From Table |^ it is also clear that the fit deteriorates quickly for Pmax > 1/6. This is a 
direct consequence of the missing kinematic constraints in the calculation (see Fig. |l|). 

To demonstrate the importance of using data that have been corrected to be for 
all particles rather than for charged ones only, we made a separate fit to ALEPH and 
DELPHI data at Mi in the massive scheme (pmin = 0.03 Mz/Q and pmax = 0.15). The 
result changes from ay^^Mi) = 0.1039±0.0009 and AiA = 0.95±0.04 GeV with xVdof = 
0.22 for all particles, to aMs(^z) = 0.1086 ± 0.0007 and AiA = 0.88 ± 0.03 GeV with 
X^/dof = 0.99 for charged particles. Thus, contrary to the case of the thrust distribution, 
the difference in as is very significant. The agreement between the values of as between 
the thrust and heavy-jet mass distributions when fitting charged particles only should be 
regarded as accidental for the reason already given. Note that is much better for all 
particles than for charged particles only. 

We have also performed fits to the experimental data using a shape function with 
exponential fall-off. For the lower cut pmin = 0.01 Mz/Q and the upper cut pmax = 0.15, 
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Figure 8: Heavy jet mass data at Mz in the decay scheme for all primary quarks together 
with the results of the fits with a shape function in the range 0.01 Mz/Q < pn < 0.15 
(solid line) and a shift in the range 0.03 Mz/Q < pn < 0.15 (dashed line). The pertur- 
bative calculation is also shown (dotted line). 



we find 

ajj^{Mz) = 0.1020 ± 0.0005 
^ , J AiA = 0.84 ± 0.03 GeV AgA^ = -0.085 ± 0.032 GeV^ 

^i-d^-^)-- 6,= 0.70 ±0.10 h^O ^^S) 

q = 0.50 ±0.77 

with x^/dof = 0.62 for 132 points. For comparison we repeated the fit in the P scheme: 

aMs{Mz) = 0.1008 ± 0.0006 
. . AiA = 0.551 ±0.012 GeV AaA^ = -0.089 ± 0.015 GeV^ , . 

f (udscbj : ^ Q ^ Q b2 = ^^^^ 

q = 10"^ ±0.11 

with x^/dof = 0.47 for 132 points. The results of these fits are in good agreement with the 
shift-based fits. This gives further support to the use of the lower cut pmin = 0.03 Mz/Q 
when fitting with a shift. 

The results of the best fits with a shape function and with a shift (with lower cut 
Pmin = 0.03 Mz/Q) in the decay scheme are shown together with the data at Mz in Fig. |[ 
The range shown is extended much above the upper cut in order to show the discrepancy 
discussed above between the calculation and the data. 

Since the precision of the resummation formula deteriorates as pu is increased, we have 
also tried to decrease the upper cut pmax when fitting with a shape function. Limiting 
the fit to include just the peak region by choosing pmax = 0.06 (the fit does not converge 
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for lower Pmax), the result is 

awsiMz) = 0.1038 ± 0.0025 
^, , AiA = 0.78 ± 0.09 GeV AaA^ = -0.057 ± 0.057 GeV^ 

D(udscb): ^^=0.45 ±0.05 b, ^ 

q = 0.33 ±0.45 

with x^/dof = 0.55 for 62 points. This shows that an analysis of the heavy-jet mass data 
restricted to the peak region is closer to the result from the thrust analysis. In addition, 
one identifies a trend for the central value of a^, which is similar to the shift-based fits: 
the lower the upper cut is, the higher gets. 

Finally, we have fitted and a shift Ai to the heavy-jet mass data with a modified 
logarithm, L = \n{l/pH — 1). The effect of this modification is very small, because of the 
limited range used in the fit. As discussed in section (^.31), this modification does not at 
all represent the magnitude of genuine non-logarithmic terms in this case. 

To summarize the independent fits to the heavy-jet mass distribution, we find smaller 
values of a^, and corresponding larger values of Ai, than those obtained in the thrust 
analysis. We understand this discrepancy by the limited applicability of the phase-space 
approximation. We have established that this primarily affects the distribution at larger 
values of ph- This is reflected for example in the fact that the best fit value of decreases 
as the fitting range shifts to the right. Owing to the limitations of the currently available 
calculations, it is not possible to extract reliably from the heavy-jet mass data. A full 
NNLO calculation can change the situation. 

3.2.2 Comparison of the heavy-jet mass with thrust-based predictions 

While the available calculation of the heavy-jet mass is not accurate far enough above 
the peak region to provide an independent measurement of a.^, the main advantage of 
DGE is that it provides a solid basis for the analysis of power corrections in the peak 
region itself. 

The purpose of this section is to test our assumptions concerning power corrections to 
the thrust and the heavy-jet mass distributions. This is done both by fitting the heavy-jet 
mass distribution with a fixed a.^, which is set according to the best fit of the thrust, and 
comparing the power corrections parameters A„, as well as by using the best fits of the 
thrust distribution to determine all the parameters and confront the calculated heavy-jet 
mass distribution with the data. 

We start by fitting a shift Ai of the perturbative heavy-jet mass distribution, for 
a fixed value of as based on the thrust fits (Table [1|). The results are summarized in 
Table together with the best fit values of Ai obtained in the thrust case for easy 
comparison. The agreement between the two observables is quite reasonable, except for 
the massive scheme, which is clearly disfavoured. We stress, however, that the heavy-jet 
mass results are quite sensitive to the fitting range and that the ^-i'^ much worse than 
in the fit where as is free. This is especially true for the data at Mz where x^/poi^t 
> 2. The sensitivity to the upper and lower cuts is summarized in Tables |ll] and |12|, 



respectively. Examining Table 11 one can get the impression that the agreement of the 



shift parameters between the two observables shown in Table 10 is misleading, since the 
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Table 10: Fits of Ai to the heavy-jet mass data with a lower cut pmin = 0.03 Mz/Q and 
an upper cut p^ax = 0.15 in different HMS using a-^{Mi) from the thrust analysis as 
input. 



HMS 


aMs(Mz) 


AiA [GeV] 


xVdof 


Points 


(AiA)thrust [GeV] 


P (udscb) 


0.1072 


0.28 ± 0.02 


1.17 


87 


0.35 ± 0.02 


E (udscb) 


0.1080 


0.38 ± 0.02 


1.32 


87 


0.42 ± 0.02 


D (udscb) 


0.1086 


0.59 ± 0.02 


1.57 


87 


0.63 ± 0.02 


M (udscb) 


0.1072 


0.85 ± 0.01 


1.15 


87 


0.35 ± 0.02 



Table 11: Fits of a shift with a fixed aMs{Mi) = 0.1086 to the heayy-jet mass data in 
the decay scheme (udscb) as a function of the lower cut = Pmin(^z) Mz/Q, with a 
fixed upper cut pmax = 0.15. 



Pmm(Mz) 


AiA [GeV] 


xVdof 


Points 


0.02 


0.70 ± 0.01 


3.16 


106 


0.03 


0.59 ± 0.02 


1.57 


87 


0.04 


0.51 ± 0.02 


0.98 


72 


0.05 


0.45 ± 0.03 


0.79 


60 



values strongly depend on the lower cut. Nevertheless, based on the shape-function fits 



presented below, the comparison in Table 11 is reliable: the shape-function-based fit 



values of Ai are consistent with those obtained with a shift, provided that the lower cut 
for the latter is pmm = 0.03 Mz/Q. 

To compare between the thrust and heavy-jet mass fits and to illustrate the correlation 
between and Ai, Fig. ^ shows the results of fitting the two parameters to the thrust 
and heavy-jet mass data in the decay scheme as 1-cr contours. The figure also shows the 
results of fitting Ai for one of the observables when fixing ag based on the other. We notice 
that the heavy-jet mass and thrust have very similar correlations between ag and Ai. 
In addition, the figure shows the results of the same exercise using NLL resummation, 
instead of DGE. Contrary to the DGE-based fits, in the NLL case there is no agreement 
between the thrust and heavy-jet mass results when is fixed. This shows that important 
contributions are missed in the NLL resummation. In addition, the magnitude of the 
non-perturbative corrections for a given is, of course, smaller in the case of DGE. 

A more stable result for the non-perturbative parameters under variations of the 
fitting range is obtained when including the peak region and fitting with a shape function. 
Fitting a shape function with an exponential fall-off in the decay scheme, where is 
fixed based on the corresponding thrust analysis and the range is 0.01 Mz/Q < pu < 0.10, 
yields 



aMs(^z) = 0.1090 fixed 
00 

g = 2 ■ 10-* ± 0.04 



AiA = 0.600 ±0.007 GeV AaA^ = 0.080 ± 0.019 GeV^ 
D (udscb): ^^^0.50 ±0.01 6,^0 ^^^^ 
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Table 12: Fits of a shift with a fixed aMs{Mz) = 0.1086 to the heavy-jet mass data 
in the decay scheme (udscb) as a function of the upper cut pmax with a fixed lower cut 
p^in = 0.03 Mz/Q. 



Pmax 


AiA [GeV] 


xVdof 


Points 


0.10 


0.64 ± 0.02 


0.95 


58 


0.11 


0.63 ± 0.02 


1.00 


60 


0.12 


0.61 ± 0.02 


1.27 


69 


0.13 


0.61 ± 0.02 


1.26 


76 


0.14 


0.60 ± 0.02 


1.47 


84 


0.15 


0.59 ± 0.02 


1.57 


87 


0.16 


0.58 ± 0.02 


1.94 


95 



Table 13: Fits of a shape function with a fixed aMs(Mz) = 0.1090 to the heavy-jet mass 
data in the decay scheme (udscb) as a function of the lower cut = Pmin(^z) Mz/Q 



cut Pmax — 


0.10 fixed. 








Pmin(Mz) 


AiA [GeV] 


A2A' [GeV2] 


xVdof 


Points 


0.01 


0.60 ± 0.01 


0.08 ± 0.02 


0.93 


99 


0.02 


0.60 ± 0.02 


-0.01 ± 0.02 


0.74 


75 


0.03 


0.57 ± 0.01 


-0.20 ± 0.07 


0.37 


58 


0.04 


0.39 ± 0.20 


-0.29 ± 0.18 


0.26 


45 



with x^/dof = 0.93 for 99 points. The reduced sensitivity of the shape-function-based fit 
to the lower and upper cuts is shown in Tables |13| and |l^, respectively. Note that fitting 
with a shape function only makes sense if the peak region is included. Therefore, the 
most relevant entries in Table |T3| are the first two lines. The agreement with the thrust 
result ( PD for Ai is very good. As for A2, although the pn fit results are not as stable 
with respect to the lower cut as those of the thrust (cf. Table 0), it seems that the pn fit 
in the decay scheme, as that of the thrust, prefers a small A2. For the parameter bi the 
agreement between the two results is also good, but this is not the case for q. It should 
also be noted that the value of q in (|5lD is at the limit of the allowed range imposed on 
the fit. 

Finally, as an ultimate test of our assumptions, we fix both and the shape function 
based on the thrust distribution best fit and examine the agreement of the calculated 
heavy-jet mass distribution with the data. This gives a x^/dof = 1.37 for 99 points in 
the range 0.01 Mz/Q < pn < 0.10. The contribution to this from each experiment 
is shown in Table ^. It is worth noting that for most of the data sets the is good 
(x^ ^1)- The only exceptions are JADE at 35 GeV, OPAL at Mz and SLD. Excluding 
these data sets the becomes x^/dof = 0.60. 

Figure ITU] shows the data at Mz together with the prediction based on the thrust 
analysis. The general agreement between the two is good. A closer look shows that the 
theory curve systematically underestimates the data up to pn — 0.06 and overestimates 
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Figure 9: The results of fitting «<, and a shift Ai to the thrust and heavy-jet mass data in 
the decay scheme using DGE or NLL resummation (l-cr contours). Also shown are the 
results of fitting Ai for one of the observables fixing based on the other (crosses). The 
ranges used are 0.05 Mz/Q < t < 0.30 and 0.03 Mz/Q < Ph < 0.15. 

Table 14: Fits of a shape function with a fixed aMs(Mz) = 0.1090 to the heavy-jet mass 
data in the decay scheme (udscb) as a function of the upper cut pmax with the lower cut 
p^in = 0.01 Mz/Q fixed. 



Pmax 


AiA [GeV] 


AsA' [GeV^] 


xVdof 


Points 


0.05 


0.63 ± 0.01 


0.08 ± 0.02 


0.58 


52 


0.10 


0.60 ± 0.01 


0.08 ± 0.02 


0.93 


99 


0.15 


0.58 ± 0.01 


0.07 ± 0.02 


1.45 


132 



it to the right of this point. In Fig. |TT| we show the prediction based on the thrust analysis 
together with the heavy-jet mass data at all energies. Here we see that the discrepancy 
between the predicted distribution and the pn data increases as the energy is decreased. 
This is natural, since the larger value of as makes the higher order corrections, which are 
unreliable for large pn, more important. 

Finally, the non-perturbative parameters determined in a shape-function-based fit in 
the P scheme are also in very good agreement with the ones found in the thrust analysis. 



cf. Eq. (^71), provided ag is fixed by the latter: 

aMs(Mz) = 0.1068 fixed 
, , AiA = 0.404 ±0.007 GeV AsA^ = 0.080 ± 0.016 GeV^ , . 
P (udscb): 5^ = 0.37 ±0.03 62 = ^^"^^ 

q = 0.46 ±0.07 



37 




Figure 10: Comparison of the heavy-jet mass data at Mz in the decay scheme with the 
theoretical prediction based on the determination of as and the shape function from the 
thrust analysis. 

with x^/dof = 0.82 for 99 points. Although the total is good, the contribution from 
the experiments at Mi is still larger than 2 per point. The result for Ai is rather stable 
with respect to variations of both the upper and lower cuts. The result for A2 turns out to 
be stable with respect to the upper cut, but, as in the decay scheme, not with respect to 
the lower cut. Fixing also the shape function to be the one found in the thrust analysis, 
the increases to x^/dof = 1.44 for 99 points. In this case the is acceptable for all 
data sets except OPAL at Mz and SLD, which contribute to 84 for 10 points and 13 
for 2 points, respectively. 

To summarize the comparison between the heavy-jet mass and the thrust, we find 
a good agreement between the two for the first moment of the shape function (or the 
shift) A 1 once Us is fixed. Although A2 is harder to determine, there is no indication for 
discrepancy between the observables even at this level. The agreement is restricted, of 
course, to the massless schemes (P, E and D). In general, the heavy-jet mass distribution 
in the peak region {pH ^ 0.1) is well described by the theoretical prediction that uses 
Us and the non-perturbative shape function from the thrust analysis as input. This 
supports our theoretical assumptions, in particular: using DGE as a basis for power 
correction analysis in the peak region (NLL resummation is insufficient), the assumption 
that the correlation between the hemisphere masses in the peak region is small (the 
simple relation between the thrust and the heavy-jet mass is based on no correlation) and, 
finally, the assumption that the non-perturbative corrections implied by the renormalon 
ambiguities dominate. In addition, we demonstrated how a careful treatment of hadron 
mass effects can open the way to quantitative comparison of hadronization effects which 
are parametrized as power corrections. 
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Figure 11: Comparison of the heavy-jet mass data in the decay scheme with the theoret- 
ical prediction based on the determination of ag and the shape function from the thrust 
analysis. Both the theoretical prediction and the data have been multiplied by a factor 
20^, with iV = ... 9 for Q = 14 ... 189 GeV. Data points with errors larger than 100 % 
are not shown. 

4 Conclusions 



Dressed gluon exponentiation is especially designed to resum large perturbative correc- 
tions in differential cross sections in the threshold region, where power-suppressed cor- 
rections are non-negligible. DGE provides the basis for parametrization of these power 
corrections. Contrary to the standard Sudakov resummation techniques, which aim at a 
fixed logarithmic accuracy (such as NLL accuracy), DGE attempts to reach power ac- 
curacy. This does not mean, of course, that the perturbative calculation itself is correct 
to power accuracy (it is ambiguous at this level) but rather that the larger, perturbative 
corrections are resummed. These corrections appear as subleading Sudakov logs that are 
enhanced by numerical coefficients, which increase factorially at large orders. 

Physically, the dominant perturbative corrections and the power corrections have the 
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same source: soft gluon emission at large angles. Although both arise from integrals over 
the running coupling, one must separate between them, resum the former and parametrize 
the latter. A systematic separation is achieved, for example, by Borel summation. The 
sum of the two contributions is physically meaningful, but each of them separately is 
ambiguous. Moreover, subleading Sudakov logs (which are above the minimal term in 
the expansion) are parametrically larger than the ambiguous power correction, and their 
functional dependence on Q and p is distinct. Having realized this, it is clear that a 
meaningful comparison with data at power accuracy cannot be achieved with a fixed 
logarithmic accuracy. DGE is necessary. 

Starting with the DGE calculation of the single-jet mass distribution in the approx- 
imation of independent emission, we expressed both the thrust and the heavy-jet mass 
distributions. Although the single-jet mass distribution is sensitive to non-inclusive con- 



tributions |39|, both the thrust and the heavy-jet mass have a reduced sensitivity: Abelian 
non- inclusive contributions appear first at N'^^^LL order while non- Abelian ones do not 
appear before NNLL (see section p.3| ). 

An immediate consequence of the resummation of running-coupling effects by DGE 
is that the Sudakov exponent is renormalization-scale invariant within a given renor- 
malization scheme. This is, of course, in contrast to NLL resummation. Fixing the 
renormalization scheme requires, in general, to fix A as well as the renormalization- 
group equation for the coupling. Here A was unambiguously fixed as A in the gluon 



bremsstrahlung scheme accounting for the singular terms in the splitting function at 
next-to-leading order, which contribute to the thrust and the heavy-jet mass at the NLL 
level. On the other hand, the renormalization-group equation was arbitrarily chosen as 
the 't Hooft scheme, in which only the first two coefficients are non-vanishing. Conse- 
quently, renormalization-scheme dependence appears (through terms that are subleading 
in /3o) at NNLL order and beyond. Let us also recall that since renormalons are resummed 
only at the exponent, renormalization-scale dependence eventually appears when match- 
ing with the NLO result. Numerically these scheme and scale dependence effects are 
small and they are sub-dominantQ with respect to other theoretical uncertainties which 
are discussed below. 

As discussed in , the DGE formula is suggestive of specific power corrections. First 
of all, the appearance of renormalon ambiguities in the exponent implies that power 
corrections also exponentiate, i.e. factorize in the Laplace space. This is compatible 
with the shape- function approach of Jl^. To a first approximation, the non-perturbative 
correction amounts to a shift of the perturbative distribution Secondly, assuming 

that renormalon-related power corrections dominate, intimate information on the shape 
function can be deduced from the singularity pattern of (pS]). The most striking property 
is the absence of a A^/(pQ)^ ambiguity in the exponent. This suggests that, if indeed 
power corrections follow the renormalon pattern, the second central moment of the shape 
function, A2, will be suppressed. 

In addition to the non-perturbative corrections that can be inferred from the renor- 
malon ambiguities, there are also finite hadron mass effects that depend on the way 
masses are treated in the measurement of the variable. We saw that differences between 
different hadron mass schemes are of the same magnitude as the non-perturbative correc- 



* Estimates have been given in section 5.6 in Q and in section 5.1 in 
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tion itself. Therefore, it is imperative that one uses the same HMS when comparing two 
variables and that this HMS is either insensitive to the masses or is defined such that the 
mass corrections for the two observables are similar. This qualifies the decay scheme as 
well as the E and the P scheme, but not the standardly used massive scheme. We found 
that, contrary to the case of average event shapes analysed in [^, the distributions anal- 
ysed in different HMS yield very similar values for a^. The differences between different 
HMS can thus be naturally recast as a redefinition of the non-perturbative parameters. 

Analysing the available thrust data, we obtained a best fit value of aMs(^z) = 0.109± 
0.001. This result is obtained both when fitting with a shift above the peak region 
and when fitting with a shape function. The consistency between the two procedures is 
reassuring. The main uncertainty in this result was found to be the effect of contributions 
suppressed by powers of t. It was estimated as ~ 4% by modifying the logarithm in the 
resummation from ln(l/t) to ln{l/t— 1). The only other important uncertainty is due to 
the HMS dependence, which was found to be ~ 2%. There is, of course, some systematic 
uncertainty from the Monte Carlo model used to transform the data between different 
schemes, which we did not investigate. Other uncertainties, which we have investigated, 
turn out to be small: ~ 0.5%, i.e. comparable to the propagated experimental error. This 
includes: using a shift or a shape function, the form of the shape function, the fitting 
range, the residual renormalization scale dependence 0, and the effect of heavy quarks. 
The use of the Monte Carlo model in quantifying the latter is not satisfactory. All in all 
our estimate for the uncertainty in as is ~ 5%. 

For the non-perturbative parameters of the shape function, our best fit values in 
the decay scheme (with all primary quarks) are AiA = 0.605 ± 0.013 GeV and A2A^ = 
0.002 ± 0.024 GeV^. The same value for Ai is obtained when fitting with shift. In other 
hadron mass schemes the value of Ai can be a factor ~ 2 lower. It is important to 
remember that the non-perturbative parameters are strongly correlated to as- In a given 
HMS, with a fixed as, the total uncertainty in Ai is less than ~ 10%. We also found a 
weak (logarithmic) energy dependence when fitting Ai separately at each energy. The 
second central moment of the shape function, A2, is correlated with both as and Ai. 
Nevertheless, it is worth noting that in the best fit in the decay scheme A2 vanishes, in 
agreement with the large /Sq renormalon pattern. This finding can regarded as a success 
of the renormalon dominance assumption. However, one should be aware that it is not 
shared by other HMS: in the P scheme A2 is significantly different from zero. Hadron 
mass effects do not have a "perturbative" origin, and therefore they certainly do not 
admit the renormalon power correction pattern. 

In the approximation of independent emission, soft gluons are radiated from a quark 
antiquark dipole. This naive phase space is useful to calculate the logarithmically en- 
hanced terms, which are associated with this configuration. However, it misses kinematic 
constraints, which are present in the exact phase space. As discussed in section p75| , this 
shows up as non-logarithmic corrections, which become important away from the two- 
jet limit. By examining our approximation versus the exact NLO calculation, we found 
that the thrust distribution is quite insensitive to such constraints, whereas the heavy-jet 
mass is. The difference between the heavy-jet mass and the thrust is well understood. 
The kinematic constraints imply correlations between the hemisphere masses: specifying 
the heavy-jet mass, as opposed to the thrust, puts a stringent constraint on the light jet 
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mass. As a consequence, the logarithmic approximation to the heavy-jet mass distribu- 
tion breaks down rather early, at pu = 1/6, and it is not as accurate as the equivalent 
calculation of the thrust distribution. Since the range of pn values where the resumma- 
tion formula is valid is rather narrow, the determination of as from this distribution is 
unreliable: independent fits yield values of as significantly smaller than the thrust and 
values of Ai correspondingly larger. We showed that the distortion of the distribution 
mainly affects larger values of ph- This is reflected for example in the fact that the best 
fit value of as decreases as the fitting range shifts to the right. Therefore, the comparison 
of non-perturbative corrections in the two-jet region between the heavy-jet mass and the 
thrust is still possible. 

Since perturbatively the correlation between the hemispheres is significant only away 
from the two-jet limit, it is natural to assume that the correlation in the peak region is 
small also on the non-perturbative level. Under this assumption, a single non-perturbative 
shape function, which is associated with the single-jet mass, incorporates the dominant 
non-perturbative corrections to both the thrust and the heavy-jet mass distributions. 

The deviation from this simple scenario was the subject of an interesting theoretical 
study and phenomenological analysis |18|. It was claimed that correlation between 



the hemispheres due to non-perturbative soft-gluon emission is essential for a consistent 
description of the heavy-jet mass distribution in the peak region^. Ref. [|I^ finds that, 
owing to radiation from one hemisphere into the other ( "non-inclusive" correction) there is 
a positive correlation between the jet masses. It should be emphasized that this dynamical 
correlation, which is introduced in [jl^ through the shape function, has a similar flattening 
impact on the thrust distribution as on the ph distribution. Whereas radiation from 
one hemisphere into the other induces a positive correlation between the hemisphere 
masses, purely kinematic considerations indicate a negative correlation, which is missed 
by the logarithmic approximation for the heavy-jet mass distribution. As discussed in 



section |2.3| , the latter seems to be the dominant feature limiting the accuracy of the 
calculation. 

We assume here that the correlation between the hemispheres in the peak region is 
small and, we confront this assumption with the data. This is done is by using the 
best fits of the thrust distribution to determine ag and the parameters of the shape 
function, and then confront the ca/cn/aiec? heavy-jet mass distribution with the data. The 
consistency of the results shows that the correlation in the peak region is indeed small. 
An alternative approach we applied was to perform fits to the heavy-jet mass data, fixing 
as by the thrust analysis. We found a good agreement between the values of the shift 
Ai in the decay and E schemes, and a reasonable agreement in the P scheme (of course, 
the values differ by much in the massive scheme). Repeating the same exercise with the 
NLL resummation instead of DGE, we found a strong discrepancy between the thrust 
and heavy-jet mass results. The agreement of Ai between the thrust and the heavy-jet 
mass persists when fitting with a shape function. For the second central moment A2 the 
results are somewhat less conclusive, since it is sensitive to the fitting range and to the 
functional form of the shape function. Nevertheless, there is no indication for discrepancy 



major differences between the analysis in |0 and ours should be noted: the resummation 
formula applied in is the standard NLL one (with an infrared cutoff), and the data for the heavy-jet 
mass are in the massive scheme. 
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between the observables even at this level. 

It is of great interest to push the accuracy of the calculations further, e.g. by matching 
with a full NNLO result when it becomes available. Some interesting questions could then 
be addressed in more detail. This includes, in particular, the consistency of the higher 
moments of the shape function with the renormalon pattern P| , the correlation between 
the hemispheres lO, ll^ , and hadron mass effects . 
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